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Abstract. We provide a thorough examination of the conditions for the 
existence of back-reaction and an "effective" acceleration (in the context of 
Buchert's averaging formalism) in regular generic spherically symmetric Lemaitre- 
Tolman-Bondi (LTB) dust models. By considering arbitrary spherical comoving 
domains, we verify rigorously the fulfillment of these conditions expressed in 
terms of suitable scalar variables that are evaluated at the domains' boundaries. 
Effective deceleration necessarily occurs in all domains in: (a) the asymptotic 
radial range of models converging to a FLRW background, (b) the asymptotic 
time range of non-vacuum hyperbolic models, (c) LTB self— similar solutions and 
(d) near a simultaneous big bang. Accelerating domains are proven to exist in 
the following scenarios: (i) central vacuum regions, (ii) central (non-vacuum) 
density voids, (iii) the intermediate radial range of models converging to a FLRW 
background, (iv) the asymptotic radial range of models converging to a Minkowski 
vacuum and (v) domains near and/or intersecting a non-simultaneous big bang. 
All these scenarios occur in hyperbolic models with negative averaged and local 
spatial curvature, though scenarios (iv) and (v) are also possible in low density 
regions of a class of elliptic models in which local spatial curvature is negative but 
its average is positive. Rough numerical estimates between —0.003 and —0.5 were 
found for the effective deceleration parameter. While the existence of accelerating 
domains cannot be ruled out in models converging to an Einstein de Sitter 
background and in domains undergoing gravitational collapse, the conditions for 
this are very restrictive. The results obtained may provide important theoretical 
clues on the effects of back— reaction and averaging in more general non— spherical 
models. 



1. Introduction. 

The dominant cosmological paradigm ( "concordance" or "A-CDM" model) provides 
an excellent fit to observational data provided an elusive "dark energy" source is 
assumed to dominate an accelerated late time cosmic evolution [1] . Alternative models 
have been proposed that challenge this paradigm under the assumption that cosmic 
dynamics follows from fully inhomogeneous solutions of Einstein's equations instead of 
perturbations on a FLRW background |lll31l31ISllSl[ZllHlini HO] • Since observed cosmic 
structure is dominated by low density regions (voids), numerous empirical models of 
inhomogeneous voids have been proposed [TTl[T2l[T3l[T4] that can fit observational data 
without assuming the existence of dark energy (see specially more recent literature 
[TSl [T51 [II]). Almost all examined void configurations are based on the Lemaitre- 
Tolman-Bondi (LTB) models [Il[inilini[lI1121l23[21[251llS], which form a well 
known class of spherically symmetric exact solutions of Einstein's equations for a 
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dust source. Although LTB models yield a very idealized description of cosmic voids, 
they have been quite useful to explore, as a rough first approximation, the effects 
of inhomogeneity in cosmic observations without the need to employ complicated 
numerical methods. As shown in |27j . the description of cosmic voids can be greatly 
enhanced by using non-spherical models. 

An important line of theoretical work among the alternative proposals to the 
concordance model is to consider the possibility that different interpretations of 
observational data could emerge from a suitable averaging procedure applicable to 
inhomogeneous sources [21 131 13 UHl IMl 1311 133 ■ While a fully covariant averaging 
procedure acting on proper tensors has already been proposed by Zalaletdinov |31i,,32] . 
the foliation dependent averaging formalism developed by Buchert, being widely used 
in the literature [5l|28l|29|, is restricted to scalars. This may be, however, sufficient to 
address the effective evolution of cosmological parameters (see [51 1331 13Z] ) • However, 
devising a well posed averaging procedure in General Relativity is still an open 
issue subjected to development and debate (see [331 1311 133 I3S1 133) ■ In Buchert's 
formalism the averaging of scalar evolution equations for inhomogeneous dust sources 
yields a variance term (the so-called "kinematic back-reaction" ) involving the squared 
fluctuations of the expansion and shear scalars. While this back-reaction term could 
effectively mimic the dynamical effect of a sort of dark energy source in the averaged 
Raychaudhuri equation, leading to an "effective" cosmic acceleration in this context, 
the interpretation of this effect in terms of actual observations is still a subject of 
debate [3 EH |38l|39] (see [40l HH |42l|43l |62] for a critical view and [4^ for a counter 
reply). 

Given their widespread use to describe cosmic voids, there is an extensive 
literature [llll71llHllini|SnilSIllS21|S31IMllSS] that considers LTB models as 
natural candidates to explore the effects and predictions of Buchert's formalism. 
The following is a quick summary of this literature. Conditions for the existence 
of effective acceleration were given by Paranjape and Singh 46J in the asymptotic 
late time evolution of hyperbolic models (as we prove in section 12, these conditions 
are incorrect). A simplified toy model was proposed by Rasanen [?71 US] for 
the interpretation of back-reaction and effective acceleration in a cosmological 
perturbative scenario and in the context of gravitational collapse. This work 
was further improved and corrected by Rasanen [49| and Paranjape and Singh 
[50] . Numeric estimations of the deceleration parameter associated with Buchert's 
formalism for particular models were presented by Bolejko and Andersson [5T] and 
Chuang et al [HI. More recently, the role of shear in estimating the magnitude of 
back-reaction was examined by Mattsson and Mattsson for central voids in hyperbolic 
models in |53j , while conditions on the sign of the back-reaction were given by Sussman 
for generic LTB models in [541 and on the effective acceleration in restricted situations 
in [Sg. 

In the present article we provide a continuation and broad generalization of |54| 
and |55| . Hence, we extend, enhance and generalize all previous literature, since we 
examine by rigorous, qualitative and numeric arguments the fulfillment of sufficient 
conditions for an effective acceleration for domains in generic models, as opposed 
to looking at excessively simplified toy models (as in [?71 HSl |SD]) or particular 
cases defined by observational constraints (as in [5T]) or suitable ansatzes (as in 
[521I53| ). We consider a wide range of scenarios in regular hyperbolic and open elliptic 
LTB models: different comoving averaging domains (central region, intermediate 
radial ranges, radial and time asymptotic ranges, near an expanding or collapsing 
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singularities). In particular, we examine the existence of effective acceleration in the 
asymptotic radial range of models converging to a vacuum Minkowskian state, which 
was not considered previously (notice that all models examined in previous articles 
[m SHI SSI [Sni El El US] are radially asymptotic to a FLRW background). However, 
we do exclude from consideration parabolic models or regions containing a symmetry 
center (because back-reaction vanishes identically [46|f52]) and models whose space 
slices have, either spherical topology ("closed" elliptic models), or lack symmetry 
centers (wormholes [20] and LTB self-similar solutions [201 ES])- 

The section by section content of the article is summarized as follows. Generic 
features of LTB models are briefly listed in section 2. In particular, we highlight the 
fact that sections of Minkowski or Schwarzschild-Kruskal spacetimes can be formally 
considered as vacuum particular cases of the models. We introduce in section 3 a 
scalar proper volume average functional acting on spherical comoving domains of 
space slices orthogonal to the 4-velocity. We also introduce an auxiliary "weighed" 
functional (the "quasi-local" average), together with the real functions associated to 
the average functionals, fluctuations and their mathematical properties. Buchert's 
dynamical equations and the sufficient conditions for the existence of a positive 
back-reaction and a negative effective deceleration parameter are given in section 
4. Since Buchert's scalar averaging is defined by proper volume integrals along the 
rest frames of the 4-velocity, we use extensively various results of previously published 
work [53 [5Hj that deals with the behavior of scalars (expansion, density and spatial 
curvature) along radial rays of space slices. We examine the sign of the back-reaction 
term in section 5, proving that it is positive in the radial asymptotic range for all 
hyperbolic models and for a wide class of elliptic models. In section 6 we provide 
qualitative guidelines describing several scenarios that should favor the existence 
of effective acceleration, identifying in section 7 specific situations in which direct 
rigorous proof of this existence (or non-existence) can be given readily without further 
assumptions (other than standard regularity). Following these guidelines we examine 
the existence of accelerating domains for central non-vacuum density voids (section 
8), in the asymptotic radial range of models converging to a section of Minkowski 
spacetime ("non-standard" and Milne sections, section 9), in low density regions in the 
intermediate transitional radial range of density clump profiles (section 10), in domains 
near (or intersecting) the curvature singularities (section 11) and in the context of the 
"spherical collapse" model (section 12). Table 4 provides a summary of these scenarios. 
In section 13 we prove that non-vacuum hyperbolic models effectively decelerate 
in their asymptotic time range and in section 14 we examine the relation between 
accelerating domains and the sign of the averaged and local spatial curvature. Section 
15 provides a summary and a final discussion of our results, together with a brief 
discussion of their theoretical implications and connections with previous literature. 

In order to keep this article as self-contained and complete as possible, we have 
included three appendices containing necessary background and support material: 
analytic solutions and regularity conditions are given in Appendix A, while rigorous 
proof of various formal results introduced in the main text are given in Appendices B 
and C. 

2. LTB dust models. 

Spherically symmetric inhomogeneous dust sources are described by the well 
known Lemaitre-Tolman-Bondi metric in a comoving frame [18] (see [19l |20] for 
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comprehensive reviews) 

r,l2 

ds^ = -di^ + -^dr^ + i?2 (d6|2 + sin^ Odcf)'^) . (1) 

where R = R{t,r), R' = dR/dr, T — J-{r) (we use geometric units G — c— 1 and r 
has length units). The field equations for ([T]) and the dust energy-momentum tensor 
T"'' — pu°'u'' with rest matter-energy density p ~ p{t, r) and u*^ ~ 5^ reduce to 

R^ (2) 

M' ^ 4tt pR'^R', (3) 

where M = M{r) and R — u'^V aR = dR/dt. The sign of — 1 determines the 
zeroes of R and thus classifies LTB models in terms of the following kinematic classes: 
J-^ = 1 (parabolic), J^^ > 1 (hyperbolic) and J-^ < 1 (elliptic). Following another 
important classification we denote the models admitting symmetry centers by "open" 
and "closed" , according to the topological equivalence class (M^ and S'^ respectively) 
of the hypersurfaces ^T[t], orthogonal to u"" and marked by arbitrary constant t. 

Besides p and R given above, other covariant objects of LTB spacetimes are the 
expansion scalar 0, the Ricci scalar "^TZ of the hypersurfaces '^T[t], the shear tensor 
fTo6 and the electric Weyl tensor E'^'^: 
^ ~ „ 2R R' . 2[(1-J"2)^i/ 

C^ab = V(,M,) - {e/'S)hab = S Eab, E^'' = U^UdC^'^" = E (5) 

where hab = UaUb + gab, Va = h^Vb, and C°'"='* is the Weyl tensor, S'^^ = h''^ - Srj'^r]'' 
with r/" = \/hP'S'^ being the unit radial vector orthogonal to u"^ and to the orbits of 
SO (3). The scalars £ and S in ([5]) are 



R R' 

R^ W 



^ = (6) 



Since the dynamics of LTB models can be fully characterized by the local covariant 
scalars {p, 9, S, £, ^TZ], the evolution of these models can be completely determined 
by a "fluid flow" description of scalar evolution equations for these scalars as in [SD] 
(see [en m [62 for the LTB case). 

We can formally include as LTB models the following two vacuum sub-cases that 
follow if the local rest mass matter-energy density vanishes: 

Vacuum-LTB models. These are hyperbolic models that follow by taking M — 
and J" > 1 in ([2|, so that we have p — from ([3]). This class of models are the 
solutions "[S2]" in (see also [15]) and describe sections of Minkowski spacetime 
in non-standard coordinates. The choice cx 1 + defines the particular case 
of the Milne Universe. We will denote by ''vacuum LTB " all Minkowski sections 
that are different from (and generalize) the Milne Universe. 

Schwarzschild— Kruskal models. If we take M = M^^^w = constant and T 
arbitrary in ([2|, then p — follows from ([s]), though a solution of this equation 
can also be given in which density is a formal Dirac delta distribution of a 
point mass associated with the "Schwarzschild mass" Mj,„hw- These models are 
representations of the Schwarzschild-Kruskal manifold in comoving coordinates 
constructed with its radial timelike geodesies (since ^ with M = Af^^h„ is the 
equation of these geodesies). The parabolic and elliptic cases are respectively 
known as Lemaitre and Novikov coordinate representation (see page 203 of [20]). 
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Notice that all covariant kinematic quantities (|4|-([6| are well defined for the 
vacuum models. These quantities correspond to a fluid flow description of sections 
of Minkowski and Schwarzschild-Kruskal spacetimes defined by a 4-velocity field 
associated by specific classes of test observers (not dust layers). 

In general, LTB models can contain regions with different kinematics (parabolic, 
hyperbolic or elliptic), and can also be matched (along a comoving boundary) 
to regions of the vacuum cases above (see [25] for a summary of these "mixed" 
configurations). The analytic solutions of ^ and regularity conditions for the models 
in general are summarized in Appendix A. 



3. Proper volume average and quasi local average. 

The hypersurfaces '^T[t\ orthogonal to a geodesic 4-velocity provide a natural and 
covariant time slicing for the LTB models. Let X{V[r]) be the set of all scalar functions 
in any compact spherical comoving domain I?[r] = x '(?[r] C ^T[i], where is the 
unit 2-sphere, t is an arbitrary constant parameter and i?[r] = {f | < f < r} is a 
semi-open set of a radial ray (spacelike geodesic) with f = marking a symmetry 
center d] We define: 

The proper volume average ( "p— average" ) is the linear functional ( )p[r] : 
X{T>[r]) — )■ M that assigns to any S £ X{'D[r]) the real number 

where dV(p) = ^ydet{hab) d'^x = J^^^R'^R' sin 6 drdOdcj) is the proper volume 
element and ..df — f^Z^ - -dr. 

The quasi— local volume average ( "q— average" ) is defined in a similar way as 
0: the Hnear functional ( ) g[r] : X (Vlr]) R that assigns to any S" e X(X'[r]) 
the real number 

(S,,M - - (S, 

Jv[r] -^dV(p) R^R'dr 

p and q functions The functionals Q and ([s]) can be used as correspondence rules 
of functions whose argument is the radial coordinate marking (for an arbitrary 
domain) the domain boundary r. We define the "p-functions" and "q-functions" 
as the real valued local functions Sp : M'^ — > M and Sq : — > M such that for all 
z>0 

Spiz) = {S)p[z], Sqiz) = {S)q[z]. (9) 

To simplify notation, and whenever there is no risk of confusion, we will use the 
symbol r as argument of these functions. The difference between the average 
functionals ( 7|) a nd (|8| and their associated functions ^ is illustrated in figure [l] 
(see also [Sllfe] ). 

f Our notation reflects the fact that spherical domains are characterized by a radial range and are 
univocally labeled by the comoving radius of the boundary f = r. This notation is easily compared 
with the standard notation: (S)[r] = {S)xi[r]- We consider domains not enclosing a symmetry center 
in sections 10 and 11. 
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It is straightforward to show directly from Q and ([8| that Sp and Sq satisfy the 
following properties (valid also for averages by replacing Sp and Sq below with {S)p 
and {S)q): 

S{r) - Sp{r) ^ --^ f S'V(p)df, S{r) - Sq{r) ^ T 5' V(,)df, (10a) 

Sp^[Spy^iS)p+iQS)p^QpSp, Sq^[Sqy={S)q + ieS)q^eqSq, (106) 

S'p = [Sp]' =^[S-Sp], S'q = [Sq]' ^^[S-Sq]^ (lOc) 



V(,)- R ^' V(,)- i?' Var)~ ^ ^ ^ 

where J-'p is the p-function associated with the scalar J- and the proper volume V(p) 
and "quasi-local" volume V(q) are given by 

V(p){r)= f dV(p)=4^/ T-^R^R'df, (11) 

Jv[r] Jo 

V(,)(r)=/ J-dV(p) =4^ / i?2i?'df= ^i?3(r). (12) 

Jv[r] Jo •J 

While functions and functionals satisfy the same derivation laws, they behave 
differently under integration (see figure 1): only the p and q functionals satisfy 
((5)[r])[r] = {S)[r], and as a consequence, they can be considered as average 
distributions complying with the variance and covariance moment definitions for 
continuous random variables EH [S5] . 




Figure 1. The difference between Sp and {S)p. The figure displays the radial 
profile of a scalar function 5(f) (solid curve) along a regular hypersurface ^T^t), 
together with its dual p-function Sp{f) (dashed curve) defined by |9|. Panels 
(a) and (b) respectively display the cases when S" < ("clump") and S' > 
("void"). The average functional ^ assigns the real number (S')[r] to the full 
domain (shaded area) marked by i9fr] = {f |0 < f < r}, whereas the function 
Sp varies along this domain. Hence, Sp and {S)p are only equal at the domain 
boundary f = r (for every domain), and so they satisfy the same differentiation 
rules (106H10C1 locally, i.e. Sp{r) = {S)p[r] and 5^(r) = (S)' [r], but b ehave 
differently when integrated along any domain. Notice, from ( |10a| and ( |10c| , that 
if 5' < in all i?[r] then S — {S)p < and the opposite situation occurs if S' > 0. 
This figure also applies for the quasi-local functions and averages. 



Back-reaction and effective acceleration in generic LTB dust models. 



7 



Applying Q to the scalars p, 9 and "^7^ in ([2|, ^ and Q we obtain the foUowing 
scaUng laws in closed analytic form (see Appendix A): 



47r 



-Pq 



M 



/Co 



i?2 



(13) 



T~La 



n 



(14) 



(15) 



i?' « 3 

where we are using (and will use henceforth) the notation 

e 3^ 

3' ^ = 

so that Hq — 65/3 and ICq = '^TZq/6 hold. By applying ([s]) to ([g]) the remaining 
covariant scalars S and £, associated with the shear and electric Weyl tensor, can be 
expressed as deviations or fluctuations of p and Q with respect to their q-duals: 

47r 

= £^-—[p~Pq]. (16) 

Given the pair of scalars {S, Sq}, it is very useful to define their relative fiuctuation 



as 



g{S) ^ Sq 



1 



^'q/ ^q _ 

ZR'/R ~ Sq{r)R^{r) 



I S'R^df, 
Jo 



(17) 



where we used the properties (10a) and (10c). Radial differentiation of Hq in (14) 
and using (10c) and (17) to eliminate radial gradients of Hq, Pq, ICq in terms of 
5^^\ S^P\ (5^'^-', yields the following useful constraint among the fiuctuations: 



2(5(«) =f],5('') + (l-r!,)(5('=). 



with f2g and I — Qq given by 
Sirpq 



mi 



Kq 



(18) 



(19) 



where we remark that VLq is also a q-scalar for being a function of the q-scalars pq 
and T-Lq (see Appendix C of [59]). The scalar fig is analogous to the FLRW Omega 
parameter, since its value depends on the kinematic class: < VLq < 1 for hyperbolic 
models, fi^ > 1 for elhptic models and fig = 1 for parabolic models. 



4. Buchert's scalar averaging formalism. 

Buchert's equations [S^ applied to LTB models (see previous literature [47l [Ml 
1151 Uni) are the proper volume average (p-average) of the 1+3 energy balance 
equation, Raychaudhuri equation and Hamiltonian constraint [§] 

{p + 3pH) ={p)-+3{p){H)=0, (20a) 

4-7r 

{ny+{nf=--{p) + 2Q[ri (206) 

{nf =^{p)-{IC)-Q[r], (20c) 

§ To simplify notation the subscript p and the domain indicator [r] wiU be henceforth omitted unless 
they are needed for clarity. However, we will keep the notation Q[r] for the back— reaction term. 
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plus the integrability condition between \20h\ and (20c 
Q[r] + 6(H) Q[r] + 2{H) {K.) + {K.) = 0, 

where Q[r\ is the kinematic "back-reaction" term given by 
Q[r]^{{'H-{n)f-{'H-n,f), 



(21) 
(22) 



and we have used ( 16 1 to eUminatc aab<^ = 6S in terms of {H—l-Lq) . Equation (20a 
simply expresses the compatibility between the averaging ([7| and the conservation 
of rest mass, but (206) and (20c) lead to an interesting re-interpretation of cosmic 
dynamics because of the presence of Q[r]. From the Raychaudhuri equation (206), we 
have the following definition for an "effective" cosmic acceleration 

A^{ny+{'Hf = --{p)+2Q[r], (23) 

together with (as is customary in the literature [Sl[25J[in]) a definition of an "effective" 
dimensionless deceleration parameter: []|| 

„ ^ ^ _ io(™' - o(2) o('") = ^""^P^ o(2) - (0A\ 



which (bearing in mind (15) and (23)) mimics the local deceleration parameter in 

FLRW dust cosmologies, with riLf'playing the kinematic role analogous to that of 
the Omega factor of a dark energy source. Therefore, considering (23) and (24), we 
shall adopt the following conventions to hold at any domain d[r\: 



Effective acceleration: 
Effective deceleration: 



^ > or q^,„. < 0, 
^ < or q^jj > 0. 



(25) 
(26) 



As a consequence, if (25) holds then Buchert's dynamical equations (20a)-(20c 



together with (23) and (24), are qualitatively analogous to the dynamical equations 
of a A-CDM model in which Q[r] is effectively mimicking the effect of a dark energy 
source. 



4-1- Sufficient conditions for hack-reaction and effective acceleration. 

Our task in this article is to examine for generic LTB models the fulfillment of condition 



(25), by looking at the sign of either A or q^^^. Unfortunately, without resorting to 
numerical methods it is practically impossible to evaluate A and q^jj, as both quantities 
involve integrals of the type ([t]) whose integrands are (in general) not given by closed 
analytic functions of (t, r). In order to proceed without having to evaluate the average 
integrals ([?]), we remark that the signs of A or q^jj in a given domain i?[r] can be 
inferred by looking at the sign of the scalars inside the (...) brackets. This can be 
achieved by using the following property that is valid for every scalar function S and 
provides sufficient conditions on the sign of {S): 

S{f) > V f e d[r] ^ {S)[r] > 0, (27) 



where we notice that the converse is not necessarily true. However, condition (27) 



is still too strong, as it requires a given sign to hold in every point of the domain 

II Since 3('H) = ^{p)/^{p) with V(p) defined by a scale factor associated with an averaging 

domain T) is often introduced [5l 1281 as a-p = V^^j , hence we can formally write: ("H) = axi/cLj) and 
■4 = ajy/axi. While this scale factor is a theoretically appealing definition, it will not be used here, 
as it is not helpful in practice to provide analytic conditions for the sign of the effective acceleration. 
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d\r\. In particular, the scalar inside the (...) brackets in (22) depends on points inside 
the domain (f) and on the boundary (r), which makes it very hard to apply (27). 
However, considering that l(H)[r] = Jipir) holds for every r, we can use the following 
identity (proved in Appendix B, see also [Ml 1551 ) 

((H(f) - mr]r)[r] = {{H{f)-H,{f)f)[r], (28) 

which holds for every domain to express the sign condition (via (27)) on Q[r] only 
in terms of the parameter f (which can take the value f — r). A sufficient condition 
for Q[r] > (which is necessary for (23)) follows by applying (27) and ( p8| to (22) 

Q(0>0 ^ Q[r]>0, with: q(r) = [H(r) - Hp(r)]2 - [H(r) - H,(r)]2, (29) 

where Q(r) is evaluated at the domain boundary f = r. 

Since we assume that p > and {H)'^ > must hold for every domain, the 
sufficient condition for the fulfillment of (25) in any domain is simply 

47r 



A > 



^> 0, 



with: A(r) = 2Q(r) 



(30) 



where (again) the converse implication does not hold, the scalars A, Q, p are evaluated 
at the domains' boundary f = r and we have used (28), (29), (30) and {'H)[r] = 7ip{r). 
Sufficient conditions for an effective acceleration can also be given by a dirnensionless 
deceleration parameter obtained by normalizing ( 30 ) with = [R/FCj^ in ( 14 ) , which 
(as T-L) is analogous to a Hubble factor. This leads (from ( 27 )) to the following sufficient 
(not necessary) condition to fulfill (25): 



q< 



qotf 



< 0, 



with: q = 



where VL^^'^ and are analogous to and rji™ in (24): 



(S) 



2 

■)(™) 



2q 

3^1 



Ho 



= no 1 + s^p^ 



Hp 

T~La 



(31) 

(32) 
(33) 



with Vlq given by (19), and we have used (18) and p — pq{\ + (5*^''^). Considering that, 
because of (27) and (28), sufficient conditions for an effective acceleration can be given 
by evaluating the scalars in ( 30 ) , ( 31 ) , ( 32 ) and ( 33 ) at the boundary of a domain i?[r] , 
it is far less complicated to verify the fulfillment of A > or q < than the fulfillment 
of > or q^,,. < from (l25|. 



5. Probing the sign of the back reaction term. 



It is useful to look first at condition (29) because it is necessary (but not sufficient) 
for (30) (which is, in turn, sufficient for (25)). Once we identify domains and models 
violating ( p9| , so that Q[r] < holds, we can state rigorously that these domains 
exhibit effective deceleration 1261). The fulfillment of condition ([29l is examined for 



domains in parabolic, hyperbolic and elliptic models separately below (see section 12 
for the mixed elliptic/hyperbolic case). We provide a summary of results in table [T] 
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5.1. Domains in parabolic models and self-similar solutions. 



The first rigorous result concerning condition (29) is that Q(r) — (and thus Q[r] = 0) 
holds identically for all domains i?[r] in parabolic LTB models or regions containing 
a symmetry center (see table [T]). The proof is trivial [52], since = 1 and 
1-Lp{r) = Hqir) hold in these domains for all r. 

Notice that Q(r) (and thus Q[r]) does not vanish in domains containing the 
"external" parabolic region of a mixed elliptic-parabolic model (see sections 12, 14 
and [M]). We also remark that (in general) ^ 1 holds in domains of the vacuum 
LTB models and of Schwarzschild-Kruskal models that are not parabolic, hence Q(r) 
is nonzero in these cases. 

LTB self-similar solutions comply with J- — J-q, where J-q > 1 or J-q < 1 for the 
hyperbolic and elliptic cases [SB] (see also Appendix C4 of [37]). From Q and (13) 
we have IC and ICq proportional to R~^. Although these solutions are not compatible 



r 



with the existence of a center worldline, and thus comoving domains like I? = S x 
are not compact, functions of the type Sp and Sq can still be constructed as in ([7| 
and ([s]) by means of improper integrals. As with the parabolic case, it is evident that 
^ = Jo 7^ 1 also implies that 'Hp(r) = 'Hq{r) must hold for all r, so that Q(r) = (and 
thus Q[r] = 0)) hold for all domains in all slices '^T[t]. However, while self-similarity 
implies JF = J-q, the converse implication is not true: as long as M (and thus pq) 
remains a free function, the condition J- = J-q defines a wider class of regular models 
(which includes self-similar solutions) for which Q[r] =0 holds in every domain 64J. 

5.2. Domains in hyperbolic and elliptic models. 

For the remaining of this article we will only consider domains in hyperbolic and 
elliptic models and regions and, unless specifically stated otherwise, we will assume 



henceforth that is not constant. Since the fulfillment of (29) in these models is 
strongly dependent on the monotonicity of the radial profiles of H and we introduce 
the notion of "Turning Value" of a scalar S (TV of S) as a point f — rtv £ "^ir] 
such that S"(rtv) =0. As a consequence, S is monotonous in 'i9[r] if there is no TV of 
S in this domain. Notice that the radial coordinate of a TV of = J-{r) is fixed for 
all slices ^T[i], while the coordinate of a TV of the time depending scalar H changes 
from slice to slice (see [55]). 



The regularity conditions (A. 13) for hyperbolic models require J- to be 



monotonous and J^' > to hold for all r (see (36) further ahead), hence these models 
admit radial profiles with either a monotonous H or with a TV ofH. In elliptic models 
the regularity conditions do not require a monotonous T (see Appendix A), hence the 
slices '^T[t] can have TV's of and/or H. Consider the following lemmas whose proof 
is given in [54] (see also [55]): 

Lemma 1. Let 'i?[r] be an arbitrary radial domain in a slice '^T[t] of a regular 
hyperbolic LTB model (including the vacuum LTB models, see section 2): 

(a) If H is monotonous in z?[r], then Q[r] > holds for all domains ^[rk] with < r, 
including an asymptotic domain covering the whole slice. 

(b) If there is a TV of H at r = rtv in a slice T[t], then values ri < rtv and r2 > rtv 
always exist such that Q[r] > holds for all i?[r] with r < ri and r > r2 

Lemma 2. Consider domains ?9[r] in arbitrary space slices ^T[t] of regular open 
elliptic LTB models: 
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(a) If Ti, and F arc monotonous in i?[r], then Q[r] < holds for aU domains 'i?[rfc] 
with Tk < r, including an asymptotic domain covering the whole slice. 

(b) If there is, either a TV oi at f = y with monotonous H, or TVs of H and T at 
f = rtv and f = y, then a value z > max(rtv, y) always exists such that Q[r] > 
holds for all ?9[r] with r > z. 

(c) If there is a TV of H at r = rtv but J" is monotonous, then values ri < min(?/, rtv) 
and r2 > max(y,rtv) always exist such that Q[r] > holds for all t9[r] with 
ri < r < r2. 

The full proof of these lemmas is provided in [51] [^and is based on using (10a) and 
( |104 to express Q(r) in ([29]) as q(r) = (f>(r)^'(r) with: 

V(p)(r) 



$(r)=H,(r)-Hp(r)= / H'{f) 

Jo ^(p) [r) 



^{r)=.2H{r)-n,{r)-n,{r)= W [f] 

Jo ^(p)\^) 



. Mr) 
Mr) 
^(P)ir) 



1 



dr, 



Mr) 



dr, 



(34) 



(35) 



Mr). 

to the monotonicity of J-p (which depends on the 



which relates the sign of Q(r 
monotonicity of J-) in the radial range '!9[r]. 

Regularity conditions (see Appendix A) imply that all vacuum LTB models 
correspond to case (a) of lemma 1 and thus Q[r] > holds in all domains of all 
their slices '^T[t]. Notice that even if there are TV's of J-' and H at f = y,rty, we 
have Q[r] < in every elliptic domain with r < min(2/,rtv)- As shown in |58j . the 
conditions for a monotonous H in the radial direction are not restrictive in non-vacuum 
hyperbolic models: it can occur at all ^T[t] or for slices ^T[t] in a wide evolution range. 
On the other hand, the conditions for a monotonous H are quite restrictive in elliptic 
models: it can only occur for some '^T[t] in a restricted evolution range (see [55]). 

Notice that taking into consideration the full time evolution range, then, in 
general, the monotonicity of H (and thus Hq) can change from one ^T[t] to the other 
in both, hyperbolic and elliptic, models (see [5S]). Therefore, the same model (whether 
elliptic or hyperbolic or mixed) can contain particular domains i9[r] in which Q(r) (and 
thus Q[r]) may have a different behavior for specific ranges of t. 



5.3. A necessary condition for a positive hack-reaction. 

It is evident from lemmas 1 and 2 and the discussion above (see also Appendix A and 
[531 [55]) that a necessary (not sufficient) condition for Q\r\ > in an arbitrary domain 
i?[r] is given by evaluating at the boundary of every domain t?[r] the condition 



J"'(r) > with F' 



3r/C, 



qO 



2T 



3rlC, 



1 „2 



2J" 



(36) 



where the subindex g in ICqo, 5q denotes evaluation at arbitrary fixed t — tQ, the 
scale factors a, F are defined by (A.3| and we have used (17), (A.4|, (A. 5) and (A. 17). 



Evidently, (36) is also a necessary condition for an effective acceleration. For domains 



in hyperbolic models (/C^o < 0) the condition (36) plus the regularity conditions (A. 13) 
imply that must be monotonically increasing in every domain dlr] (and in every slice 
^T[t] because of (A. 5)). For domains in elliptic models (/C^o > 0), we have ^(0) = 1 



^ The proof of the case with TV's of and H in 54 (proposition 7) has an error: since we need 
to consider domains for which r > max(?/, rtv) holds, then the possibility that Q(r) > in the range 
r < min{j/, rtv) can be ruled out. 
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Domains in parabolic models: J- = \ {ICq = 0) 


Turning values 


Domain restrictions on Q[r] 


Comments 


None 


Q[r] = holds for all i9[r] 


See I11IS2] 


Domains in hyperbolic models: J- > 1 {ICq < 0) 


Turning values 


Domain restrictions on Q[r] 


Comments 


None 


Q[r] > holds for all d[r] 


Lemma 1(a) 


H 


Q[r] > holds for d[r] with 
< r < T"! and r > r2 with ri < r2 


Lemma 1(b) 


Domains in open elliptic models: J- '' < 1 {ICq > 0) 


Turning values 


Domain restrictions on Q[r] 


Comments 


None 


Q[r]<Q holds for all i9[r] 


Lemma 2(a) 


n 


Q[r] > holds for i)[r] 
with ri < r < r2 


Lemma 2(c) 




Q[r]>Q holds for 79[r] 
with r > y with y > 


Lemma 2(b) 


F and % 


as above 


Lemma 2(b) 



Table 1. Summary of sufficient conditions for tlie sign of back— reaction 
in LTB models. The conditions are expressed (second column) in terms of the 
sign of Q[r] in domains characterized by the radial range i?[r] = {f\0 < f < r} 
for each class of models, according to the existence of turning values discussed 
in section 5 (first column). The third column lists the lemmas dealing with each 
case. The domains correspond to arbitrary regular slices ^T[t] not intersecting 
a singularity (sec figures 8 and 9a). The parameters ri, r2 and y depend on the 
profiles of and H for each case and are given in the proofs of lemmas 1 and 2 in 
[54 1 . See section 12 for a discussion of back— reaction in mixed elliptic/hyperbolic 
and elliptic parabolic configurations. Back— reaction for closed elliptic models is 
discussed in |54| . 



and J^' < near r = 0, hence a TV of J' (case (b) of lemma 2) implies that J^' > 
holds for r > rtv and — 1 as r — ^ oo. Hence, Q[r] > necessarily holds in domains 



in which (36 1 holds and spatial curvature is rapidly decreasing (see section 13). 



6. Probing the existence of effective acceleration. 

In domains in which we have proven that Q(r) > or Cl'-'^^ > hold, the conditions for 
an effective acceleration (30) or (31) become a comparison of two positive quantities 
evaluated at the boundary of any domain i?[r]. However, it is not (in general) a simple 
task to verify in these cases if the positive or negative contribution dominates in this 
comparison. 



6.1. Qualitative guidelines. 



Assuming domains such that £7'^'^^ > holds (positive back-reaction) in (32), then the 
following two conditions (which arc not mutually exclusive) yield the scenarios with 
the best possibilities for (31) to occur: 



is as small as possible. This condition occurs in the following scenarios: 
(i) Strict density vacuum: p pq — flq — Cl'^''"^ — 0. See section 7.2. 
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Schwarzschild density vacuum: p — Q but pq > Q 
have fJ'™) = but ^^„ > 0. See section 7.2. 



(ii) 

(iii) Very low (but nonzero) density: ICq < and pq <^ \IC, 



1 + S'-P'^ = 0, hence we 



fig <C 1 (see 



(19 1, (A. 5) and (A.6|). Domains in hyperbolic or asymptotically hyperbolic 
models in the following conditions: 

t finite: central void region (see section 8). (37a) 
t finite: asymptotic radial range (see sections 9 and 10). (376) 
i — >■ oo asymptotic time range (see section 13). (37c) 



(iv) 



oo 

l + d^P'> 0. 

or regions. From (|17|) applied to 



This scenario is compatible with hyperbolic and elliptic models 
S — p, it occurs if p ^ under the 



following conditions of very small and very large density: 

{p,Pq} — radial asymptotic range (see section 9). (38a) 

{p,Pq} — > oo near curvature singularities (see sections 11 and 12). (386) 

is as large as possible. This condition implies the largest possible growth of 
back-reaction, which couples to density and spatial curvature [5) J9 through 
their radial gradients. The necessary condition for Q[r] > in (36) shows this 
coupling with the gradients of spatial curvature through the relative fluctuation 
S^^f which relates to the gradient IC'q through (17 1. The relation between the 
growth of back-reaction and the existence of large radial gradients of density 
and/or spatial curvature can also be appreciated from the fact that Q{r) in (29) 
involves the fluctuations H — Hp and H — Hq, which relate to the gradients Hp 



and H'q by means of (10c). In terms of these gradients we obtain (see [55]): 



= 2 



'^p/'Hp 



"-q 



- 2 



(39) 



with T-L'p and H'^ related to gradients of density and curvature by means of the 



constraints (14) and (20c). However, the coupling between gradients of density, 
spatial curvature and expansion can be complicated, and thus there may be several 
scenarios in which these gradients could "cooperate" so that US'^^ increases and 
decreases. A favorable scenario may be a low density region (small f2'^™^) in 
which fluctuations of spatial curvature are larger than density fluctuations (see 
sections 8 and 10), so that the growth of Q(r) is mostly driven by the gradients 
of spatial curvature in that region. As commented by O [29l Ell |48]), large 
fluctuations of spatial curvature occur when the functional form of (/C) = /Cp is 



different from a "FLRW-like" dependence (oc where 



7^ - 
'"D 



V(p))- 



It is important to remark that the scenarios we have described above simply provide 
the guidelines that may favor the existence of an effective acceleration in terms of 
sufficient (not necessary) conditions that are evaluated in the domain boundary, and 
thus we cannot rule out such acceleration arising under different conditions. Notice 
that even when these guidelines are satisfied the actual existence of this acceleration 
must be verified by looking carefully at the behavior of the scalars in (30) or (31) in 
the appropriate domains. 
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7. Scenarios allowing for a direct proof of effective deceleration and 
acceleration. 



Consider the following limits along radial rays of a fully regular and complete slice 
^T[t] with arbitrary finite t: 



lim S{r) = Sc = 5(0) 
lim S{r) = 5*^ 



<^ lim Sp{r) — lim Sn{r) = Sc, 
lim Sa(r) — S , 



(40 a) 
(406) 



lim Sp{r) 

where Sc and are real constants. If 5 = S{t, r) and Sc and S^ are nonzero, then 
the latter are (in general) different constants at different slices ■^T[i] (i.e. we have: 
Sc — Sc{t) and S^ = S^{t)). The subindex c will henceforth denote evaluation at 
r = 0, while ^ will denote the asymptotic limit as r — > c». By applying (40a I and 



(40 6 1 to the definitions (30) and (31) we obtain direct proof of existence of effective 



deceleration or acceleration in specific situations. 



7.1. Effective deceleration. 



Non vacuum center and central region of elliptic models. If pc > then (14 1 
implies He ^ (with K-c < Snpc/S). From ^ we have ^(''^(O) = ^^^^0) ^ 
5(«)(0) = and (40a) implies HpiO) = H,(0) = He, hence we have an effective 



deceleration in the center worldline of every regular model with a non-vacuum 
center: 



Q(0) = 0, 



A(0) = 



47r 



-Pc < 0, 



q(o) 



> 0. 



(41) 



3 , -xw 2 

Since Q[r] < holds for domains sufficiently close to the center in all elliptic 



models (all cases of lemma 2), then (23) and (26) imply effective deceleration in 
these domains. In domains of hyperbolic models Q[r] > holds for r > near 
the center (in the two cases of lemma 1), hence it is not immediately evident if 
(231 and (25) hold for r > 0. We examine this possibility in sections 7 and 9. 
Radial asymptotic range in models converging to FLRW. LTB models can 
be radially asymptotic to spatially flat (Einstein de Sitter) or negatively curved 
FLRW models, but not to positively curved FLRW models (see details in [57]). 
LTB Models wi th th is asymptotic behavior are characterized by > 0, and thus 
(14 1, (17| and \40b\ imply for expanding models > and S^^'> ^ S^^'> = 0, 
wmle n 



and Hq as r ^ oo (see |S7j)- Hence 



lim q(r) = 

r— ^oo 

where and fl 



An 

lim A(r) = — -p^ 

•->-oo j 



< 0, lim q = q — 



>o, 



(42) 



"background 
Evidently, 
range of these 



correspond to the density and Omega factor of the 
FLRW state reached asymptotically along the radial rays [57) . 



( 42 ) implies an effective deceleration ( 26 ) in the radial asymptotic 
models. 



The same situation occurs in the radial range of 
"Swiss Cheese" configurations in which the FLRW background is not reached 
asymptotically but is forced by a matching with an LTB region in a comoving 
boundary marked by finite r. We remark that all models examined in previous 
literature [47l |48l |49l [50l [511 [52l [53] are radially convergent to FLRW cosmologies, 
hence there is effective deceleration in domains extending into their asymptotic 
radial range. 
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(a) (b) 




Figure 2. LTB v acu um region. The figure displays the ratio Hp/'H (panel (a)) 
and q defined by ( |31| l (panel (b)) for an extended central vacuum LTB region. 
The region can be matched to an expanding non-vacuum hyperbolic model at 
any fixed r. The curves (from top to bottom in (a) and bottom to top in (b)) 
correspond to t = 0, 5, 10, 20, 100, with t = to = marking the fiducial slice where 
we define the initial value functions Pqo = and Kqo = 1/(1 + r^/"*). Notice how 
both plotted quantities tend to zero for large t and tend to an asymptotic value 
for large r. Panel (a) shows that the averaged Hubble factor (whose magnitude 
is Hp) can be up to 30% larger than the local Hubble factor. 



\ 

\m(r) 
\ 

\ 

\ 

2Q(rjV^ 


N(ri>>^ 


k(r) 





LTB interior Schwarzschild exterior 



Figure 3. LTB inner region matched to a Schwarzschild exterior. The 

figure display s qu alitative plots of the functions 2Q(r),m(r) = 47rp/3 and A(r), 
appearing in (ISOjl, at an arbitrary slice ^Tlt] of a configuration made of a region 
of an LTB non-vacuum model (shaded area) complying with Q(r) > matched 
with a Schwarzschild exterior at r = r;,. The local density complies with p = for 
r > ri,, hence for this radial range we have A = 2 Q > 0, though A may be negative 
in most of the inner region where p > 0. 
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7.2. Effective acceleration in vacuum regions imposed by matchings. 

Vacuum central region. This is the strict vacuum scenario of case (i) in the 
guidelines above: a vacuum LTB model extending (see section 2) between r = 
and a given finite r = rf,, can be smoothly matched with a section of a non- 
vacuum hyperbolic model extending for r > rj, (see |25|). Matching and regularity 
conditions imply that p{rb) = must hold and p must be continuous [25] . which 
implies that density must have a radial void profile {i.e. p' > 0) for r > rf, 
for all the time evolution. However, as shown in [S^, a void profile for all the 
time evolution is only compatible with absence of shell crossings in non-vacuum 
hyperbolic models if they have a simultaneous big bang (t[_^ = 0) , and this feature 
necessarily implies (from [ST]) radial asymptotic convergence to a FLRW model. 
Since p{r) = and (from lemma 1) we have Q{r) > in the vacuum region, 
then A(r) = 2Q(r) > and q(r) = —Cl'-'^'>{r) < hold in this region, and thus, 
( [25| holds by vituc of ( [SO] ). While A(r) may remain positive near rt,, the radial 
convergence to FLRW implies that A < holds in the asymptotic radial range (as 
in (42 )). We plot in figure[2]the deceleration parameter q and the ratio 'Hp/'H for a 



central vacuum region (the matching with a non-vacuum region is not displayed) . 
As shown by the figure, both quotients tend asymptotically to constants whose 
maximal values are q ss —0.03 and Hp/H ~ 1.35. 

Schwarzschild exterior. This is the Schwarzschild vacuum scenario of case (ii) in 
the guidelines above: a spherical comoving region of any LTB model enclosing 
a symmetry center (the "interior") can be smoothly matched at r = rf, with 
a section of a Schwarzschild-Kruskal manifold (the "exterior") (see page 332 of 
[5D]). Notice that pq = holds only for vacuum LTB models (M = 0), but not 
for Schwarzschild-Kruskal models: (47r/3)p(j = M^^i,„/R^. In particular, we can 
obtain this configuration by choosing the local density so that p{rb) — 0, with 
p{r) = for r > rb- Hence we have S'^^^rb) = —1, and p — together with 
= -1 for r > rt, while (ITsl) implies that (47r/3)pq(rb) = 2M{rb)/ R^{rb) > 



holds. However, the smoothness of the match [25] does not require the interior 
density vanishing at r^, so that p and S^''^ can have a step discontinuity at the 
interface, though pq, )Cq and Hq must remain continuous at rb. If Q(r) > 



holds in the whole inner region (see lemmas 1 and 2 and condition (36)), then 
depending on the density, A may not be positive for all r < r^. However, since we 
have p = in the Schwarzschild exterior, then A(r) = 2Q(r) > will necessarily 



hold at least for parts of the range r > rb (and thus (|25|) holds by vitue of (30)). 
Evidently, there will be an effective acceleration in radial ranges of the exterior 
where Q(r) > holds, and thus (following lemma 1 and lemma 2), these ranges 
depend on whether the exterior is hyperbolic or elliptic (we omit the case of a 
parabolic interior because Q(r) = 0). Figure [3] provides a qualitative depiction of 
A(r) for an LTB region with Q(r) > matched to a Schwarzschild exterior. 

It is evident that an inner vacuum LTB region or an outer Schwarzschild exterior are 
rather artificial configurations. Probing the existence of accelerating domains in less 
artificial cases is obviously more difficult. 
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Figure 4. Expanding non— va cuu m void. Pannel (a) displays the functions 
AC), f2('") and -q, appearing in for varying r in tlie initial slice t = = 

for a non-vacuum central void in a hyperbolic model. The functions were obtained 
from a numeric integration of the implicit solution l |A.8[ l for initial value functions 
that yield an evolution free from shell crossings (see Appendix A). Notice that 
as r grows back-reaction tends asymptotically to a small value {tl^^'l ft: 0.0006), 
while local density (Si'"'') grows much faster. Hence, q is negative only in a 
small radial range. Panel (b) displays q in for the values (bottom to top) 
t = — f , 0, 1, 2, 3, 4. Notice that q — !> as t grows. 



8. Under dense void regions of hyperbolic models. 



While ( 41 ) implies a positive q — at the center worldline if pc > it is possible 
to find conditions that facilitate a negative q for finite r > in non-vacuum regions 
with low density. Since regularity and radial void density profiles are only compatible 
for hyperbolic models 58J, we consider an expanding hyperbolic central void region 



having a low (but nonzero) central density, so that f2g(0) <C 1 (case (37a)) 

Following the qualitative guidelines outlined in section 6, a 
conditions that may favor the joint fulfillment of a small US™"'^ 
a range < ri < r < r2 are: 



plausible set of 
and a large 17''^^ in 



1 - 



e < 1, 



1, 



< and (5^^' < 0. 



(43 a) 
(436) 



and clearly correspond to the low density scenario (37a), together with a growth of 



ri'''^^ that is mostly driven by the variation (gradients) of spatial curvature (|(5*^'^^|). 
On the other hand, (43 6 1 follows from the regularity conditions, since a regular void 
58 that (5^''^ > and (5*^'^-' < hold, and is also compatible with the 
(jCC) ^ Q (-vvhich implies %' < 0, leading from (34) to 



profile re quires 
fact that (43a) implies 26^'^'> 
a negative 1 — V.p/'Hq). 



Under the assumptions (43a) and (436) condition (31 1 takes in this radial range 



+ The possibility of having pc = only in the central world line is equivalent to reducing the vacuum 
LTB region of the previous section to this single worldline. For non-vacuum models this situation 
can be considered as the limit pc —> of very low (but nonzero) central density. While He = 
can arise mathematically if either < ICc = (87r/3)pc or Kc = pc = occur, both possibilities are 
incompatible with regularity conditions and will not be considered. 
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the approximated form 



Hp 
T~L„ 



Hp 



+ o <0, 



(44) 



which may hold if e/2 is effectively smaller than the other term. Evidently, the chances 
for (44) to hold improve if the non-vacuum void is "deeper" (smallest possible e ^ 1) 
and the larger we have the ratio between the gradients o f den sity a nd c urvature 
(|4''V4'^^I ~ 0)- The "wider" we can set up the void (( [ii^ and ( |436| holding 
for the largest possible radial range), the larger radial range may exhibit q < 0. In 
the limit pc => fig(O) — >■ 0, the void region that fulfills (44 1 extends all the 
way to r = (where q — >■ 0). Since it is very difficult to obtain analytic constraints 
on 1 — Hp/'Hq, we provide in figure |4] a numerical example of a central void region in 
which q < holds in an intermediate region close to the center. 

Models of non-vacuum voids were examined by Mattsson and Mattsson in [53] , 
who explored the relation between the growth of back-reaction, the average of the 
shear scalar and the smoothness of the gradients of density. As argued in this reference, 
the use of discontinuous step function density radial profiles made by joining regions 
with constant density (for example glueing FLRW regions as in [47l |48]) yields a 
gross error in the estimation of the magnitude of the back-reaction in comparison 
with smooth profiles, since shear is zero in each constant density region. The 
connection between back-reaction and shear follows from the fact that the shear scalar 
CTafecr"^ = 6I]2 = 6(-H/(«))2 (from ^ and ^) is related to H'g through ^ and 
Hence, shear may favor an effective acceleration in a given radial range in low density 
non-vacuum voids if it is associated with large gradients of spatial curvature. 



9. Effective acceleration in the radial asymptotic range. 

Since there is no effective acceleration in domains extending into the radial asymptotic 
range of LTB models that converge to non-vacuum FLRW spacetimes (section 7.1), 
we consider in this section only models converging in their radial asymptotic range to 
asymptotic states that are not FLRW [57] : 

Models radially asymptotic to non— Milne sections of Minkovirski (vacuum LTB). 

These hyperbolic and elliptic models are characterized by the asymptotic limits 



{p^,JC. H^} = [STj. Since we have Hq ^ and Hp — >■ 0, we cannot use (17) 



and (406) to obtain the asymptotic values of (5*^''^ and 1 — Hp/Hq in this limit 
The asymptotic limit of flq depends then on the asymptotic behavior of the ratio 



1^ I ( 0, ^q^ 1, Matter dominated (MD) 
x= ^ I x^, nq^ilq^, Generic (G) (45) 

3 Pi [ oo, ^ 0, Vacuum dominated (VD) 

where in the G case ilq^ > 1 if /Cg > and < ilq^ < 1 if /C, < 0. Notice 
that for all elliptic models we must have ICq < (87r/3)pq, hence x < 2 holds and 
all these models are either MD (x^ < 2) or G {x^ ~ 2). On the other hand, 
hyperbolic models can be MD, G or VD. 

Models radially asymptotic to Milne. These models must be hyperbolic and 
VD, corresponding to [57] pq ^ — and JCq JC^ < 0, and thus 
•H^ = |/C^|i/2 and ^ 0. 
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While some results can be obtained without making specific assumptions on the radial 
asymptotic convergence of the involved scalars, it is useful (and necessary) in most 
cases to consider these assumptions to prove the existence and an estimated numerical 
value of an effective acceleration in this range. Following reference [57j, the radial 
asymptotic forms of all scalars can be obtained from the asymptotic behavior of the 
initial value functions pqQ and TCqo (evaluated at an arbitrary t ^ to). In particular, we 
will consider convergence to the following power law trial functions (the same results 
were obtained with other trial functions in f57J ) : 

S^"''--^, (0<a<3), 

< /3 < 2, (hyperbolic models) 



Pq ^ Pqo ^ mo r' 



(46) 



qO 



kor 



(3_ 
3' 



/3>2, (elliptic models) ,(47) 
(3 — 0, (asymptotically Milne) 

where the symbol ^ indicates uniform convergence (see |57j for the precise definition) , 
mo, ko are nonzero real constants with inverse square length units, and the ranges of 
the parameters a and /3 follow from compliance with the regularity conditions (A. 13 1 
and (A. 14). In particular, the range of /3 in elliptic models follows from T with ICqo > 
in (A.4|. Models asymptotic to self-similar LTB solutions are the G models complying 
with a = (3 — 2, so that the limit (501 holds, while models converging to Schwarzschild 
are characterized by a = 3 (so that M — >• M^^^w, see [S7j). 

Since the free parameters commonly used in the literature are the functions M 
and T (commonly given as "J" = \/l + 2i?" ) in ^ and it is useful to provide the 
asymptotic form of these parameters under the assumptions (46)-(47): 

M^mor^-", [l-fcor^-^]'^', (48) 

while the bang time follows by applying (46) and (47l to (A.15I, leading to 
tbb ~ l/\ko\^^^ for models converging to Milne, while for models converging to a 
vacuum LTB state we have 

t,, ^ to - mo|A:or3/V3^/2-"Z (|fcol r"-^/mo) , 



(49) 



where Z stands for or Zfi defined by (|A.10|) and (|A.11|) for hyperbolic or elliptic 
models, 



so that —oo as r — )■ 00 [57j . 

As opposed to models converging to FLRW, 



finding the asymptotic behavior 



of q now requires further examination, as we cannot use the limits (406) and (41)- 
(42) to infer the asymptotic limit of the ratio 1 — Hp/Hq. Whether we consider the 
assumptions (46)-(47l or not, we remark that the difference between Hg and Hp is 
the presence of the function in their definitions ([7| and ([s]) (see also (34)), hence 
the limit of this ratio depends on the asymptotic behavior of this function. As we 
prove in Appendix C, the following two cases arise 



lim J- 



lim jF = oo 



lim 1 



lim 1 



Hp 



0. 



(50) 
(51) 



where the form of ^ > depends on the asymptotic convergence of pq, JCq, Hq (see 
(C.9)). Notice that both limits (50) and (51) are compatible with the necessary 
condition for Q[r] > in (36). We can now infer the existence of an effective 
acceleration associated with q — )■ q < according to the following cases: 
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Figure 5. Value of the effective acceleration in the radial asymptotic 
limit of G and VD hyperbolic models. The asymptotic limit of the 
function q in l |57| for hyperbolic G models (case (C)) is displayed as a function 
of the parameter 7 given by | |57| l for values (bottom to top) (loo = f^<joo = 
0,0.011,0.022,0.033,0.044,0.055,0.066,0.077,0.088,0.1. The thick curve is the 
case = 0. Since l |57[ l reduces to | |54[ l, this case corresponds to VD models 

(case (A)). Notice that a very small value flq^ < 0.09 is needed for having q < 
in the G models, while for the VD models q reaches a maximal value q ~ —0.035 
for 7 = /3 Ki 1.25. 



Case (A): there is efTective acceleration. The only case in which a conclusive 
proof of existence of an effective acceleration can be given without making 
assumptions on asymptotic convergence are the class of LTB models complying 
with: 



ICq < (hyperbohc), 



J" -> 00 [eq. (51)], 



(VD). 



Independently of the assumptions (|46|)-(|47|), condition (|3l|) takes the limit 
lim q = q 



<0, 



(52) 



(53) 



where J^-* is the asymptotic limit of 6'^'^\ which must be negative since 5'^'^-' < 
must hold asymptotically because \ICq\ 



and thus 



< 0. This case 



corresponds to the asymptotic forms (46) and (47) with 3 > a > 2 > /3 (models 
converging to Schwarzschild are the particular case a = 3). Considering ( |52[ ), we 
have n 



the limit ( 53 ) 



/3/2 from (C.9), so that (31) and (46)-(47) imply the following form for 
/32(2-/3)(10-^) 



288 



(54) 



whose plot in figure [5] (the thick curve) reveals that it has a maximum of 
q « -0.035 for /3 « 1.25. 

Case (B): there is no effective acceleration. Conclusive proof that q -H- q^ > 
holds can be readily given for LTB models that are: 

(i) asymptotic to FLRW dust cosmologies (see section 7.1). 

(ii) asymptotic to vacuum LTB that are elliptic and violate the necessary 



condition (|36) 



2), and thus 

< 1, so that J^' < and Q < hold asymptotically (case (a) of 



Hence, (from (17)) Kq decays a,s R {P 
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lemma 2). Models asymptotic to elliptic self-similar solutions are contained 
in this class. 



(iii) asymptotic to vacuum LTB that are elliptic and comply with ( 36 ) , so that 
Kq decays faster than (/3 > 2), but pq decays slower than {a < 3) 
from (jls])). Hence 1 — ILjHq — > and Vlq > 1 hold, but we have 
5^P^ -^^^ > -1 (from mi). 

(iv) asymptotic to vacuum LTB that are hyperbolic with T > 1 (case 
( [50| ) and that are MD and G, since 1 - "Hp/H, with Qq^ = 1 (MD) 
and < fig^ < 1 (G), and S'^''^ — ?> S^^^ > — 1 because pq decays slower than 
ICq (since a; — > or x — > in (45)). Models asymptotic to hyperbolic 
self-similar solutions are contained in this class. 

Case (C): there may be effective acceleration (Limit of q is the difference of 
two positive quantites) in models complying with 



/Cq < (hyperbohc). 



Qq 



< r!„ < 1, 



(55) 



(G), 
(MD), 

which arc the MD and G models of case (b) above, and lead to the exact 
asymptotic limit: 



(56) 



where S^'^'> and 5'^''-' must be nesative because p' < and K.' < must hold 



asymptotically to reach pq ^ and ICq 0. The sign of the limit of q in (56) 
depends on the models being MD or G: 

G models: we have a = (3 = -f < 2, with (56) taking the form 
7^(2 



7)(10-7) 
288 2 



1 - 



2mo 



2mo 



< 1, 



(57) 

which contains the VD models as the limiting case ilq^ 0. The graph 
of q^ as a function of 7 for various values of flg^ is displayed by figure [5j 
revealing that q^ < holds only for Qq^ < 0.09 and taking the largest 
values (and ranges of 7) in the limit V,q^ — >• (case (A) above). 
MD models: equation ([56]) takes the form: 
, _ a2(/3-2)(/3-2a + 10) 1 

~ 18(/3-a + 4)2 +2V" 6. 

whose sign can be found by plotting the right hand side of this equation as a 
function of a, /3. Such plot (not displayed) reveals that there is no effective 
acceleration for these MD models, since q^ is negative only in a range of 
values for which a > {3. 

Case (D): there may be eflfective acceleration if q < as q — > 0. This occurs 
in LTB models that are: 

0. 



a< (3 <2. 



(58) 



(i) asymptotic to the Milne Universe, since 1 — 'Hp/'Hq — > but fiq 

(ii) elliptic asymptotic to Schwarzschild. These models comply with (36 1, with 



ICq decaying faster than (/3 > 2), so that there is a TV of F and 

— >■ 1, Q > hold asymptotically (case (b) of lemma 2), but pq decays 
as (a = 3), so that 1 + 5^^^ — >■ holds asymptotically. Notice that, from 
(13 1, (46) and (A.4|, we have in this limit M — > TOoi?Q = const, which can 



be identified with an asymptotic Schwarzschild mass ilC,hw [57] . 
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(iii) hyperbolic with F — > > 1 (case ( |50[ )) and that are VD, since we have 
1 - Tip/nq but also flq-^O holds. 



9.1. A closer look at models of Case (D). 

We examine below the possibility that q < holds as q — > for models asymptotic to 
Milne (case (D)(i)) and elliptic models converging to Schwarzschild (the case (D)(ii)). 

Models converging to Milne. These models are characterized by asymptotic 
forms in (46) and (47): Pq ^ mor^°' and \ICq\ ~ fcg > (/3 = 2) up to the 



first leading term. We consider their asymptotic forms up to a second order term: 



Pq ^ mgr + mir 



a < ai, 



\)Cq\ ^ ko + kir 



-Pi 



(59) 



where ki and /3i are positive. Following the guidelines outlined in section 6.1, the 
best conditions for q < (lowest possible near constant density and significant 



gradients of /C,), we assume that a> Pi. Hence, we get from (14), (18), (191 and 



(481 



T-Ln 



2ko 



1 



k^ 
2ko 



3fco 



mo _Q 
mo 

F « \/kor, 



(60a) 



(606) 



where we are using the symbol w instead of ~ because the convergence may not 
be uniform. In order to compute the asymptotic form of Hp, we follow the proof 
of the limit (51 1 in Appendix E to evaluate the following integrals from ([t]) and 

(Inl: 



2n 



- V(rt(y) + — r^ V^p^iy) = V^p^iy) - 



An 



nv; 



(9) 



/fco 

df Ip{y) + 27rr^ 



where: Ip{y) = in —dr — 2nr 

Jo 

From ([7]) and considering r ^ y we have 



fci(l-/3i/3) 
2fco(l-/3i/2) 

, fci(l-/3i/3) 



2fco(l-/3i/2)' 



-/3i/2 



fci(l-/3i/3) ^ 
2A:o(l-/3i/2)' 



Hp 

T~La 



kiPi r 



2fco(l-/3i/2)' 



which leads to 

„ fcl/3i(3-/?i) _2^^ 

6fco(2-/?i)2'' 



2ko 



1 



a 
3 



mi{ai — a) 



Too 



(61) 
(62) 
(63) 

(64) 
(65) 



Hence, condition (31) is fulfilled in this case if /3i < 3 and a = 3 hold (since 
ai > a). It is also fulfilled for a < 3 as long as /3i < 3 and a > 2/3i hold, which 
implies that /3i < 3/2 must also hold (since a < 3). 

Elliptic models converging to Schwarzschild. In this case S^^^^ 1 implies 

Pq cx r~^, we consider the following asymptotic convergence forms up to second 
order 



Pq ^ TOor 



TOir 



ai > 3, 



kor 1^, /3 > ai, 



(66) 
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-3/2 



leading to 
Ti.q « -\/2mo r 

2^(«) 



1 - 



4mo 



J"! 



3) 



,3-ai 



3too 



3' 



Considering that % 



kg 

2mn 



(67a) 
(676) 
(67c) 



^^(1 + S^^^) and using (181, we construct the same type 



of integrals as in (61 1-(63) to compute Hp, leading to 

15/4) 



mor 



-3/2 



1 



mi(ai 

2mo(9/2 - ai) 

so that using (67a)-(67c) and (31), and bearing in mind that /3 > ai, we get 

3)2 



(68) 



mi{ai 



3— ai 



ai) 



(69) 



^ 6mo(9/2 

which is negative if ai < 9/2. If we assume that 3 < /? < ai we obtain 
q cx ~r^^^ > 0. In either case, the value of q is much smaller than in the 
hyperbolic models in which < 0, as can be appreciated by comparing the 
curves of q in figures 7a and 7b. It is important to remark that the fact that /3 > 3 
implies that ^q'^'' < — 1 must hold asymptotically (from (36)), hence as we show 
in section 13 (see figure lib), local spatial curvature JC —^^/6 is asymptotically 
negative for these elliptic models. 



10. Effective acceleration in the radial intermediate range of density 
clumps. 

If we have a clump density profile, either in domains of hyperbolic or elliptic models 
or mixed elliptic/hyperbolic configurations, we can think of an "intermediate" radial 
range associated with the transition between "near center" and asymptotic conditions. 
Evidently, the radial gradients of the basic scalars p, /C, H take their largest value in 
this region, and thus we can apply the guidelines specified in section 6.1 that favor 



accelerating domains like (43a) and ( |43&[ ). Following the same assumptions described 
in section 8 for void profiles, the combination of small Jl^™) and large Cl'^'^^ occurs in this 
range if there are large density and/or curvature gradients, but now the low density 
follows as a sharp density decay from a large central value, remaining sufficiently 
small in the full range r > r2. Hence, the possibilities for mid-range accelerating 
domains strongly depend on asymptotic convergence of the models. We examine 
these possibilities for models converging to FLRW and vacuum LTB states. 

Models converging radially to FLRW. While q eventually becomes positive for 



sufficiently large r (from (42)), it can be negative in a low density region in the 
intermediate radial range. The conditions for a small Jl^™) are sensitive to both 
1 + S^P^ and Qq. We have the following cases (see table 3 and [57]): 
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Figure 6. Effective acceleration in the intermediate radial range: 
models converging to FLRW. The figures depict the deceleration parameter 
q defined in ( |31| (thick solid curve), together with the Omega factors for back- 
reaction, and local density, O'™' (dashed and dotted curves), for a fiducial 
initial slice t = to = oi sl mixed elliptic/hyperbolic models that converge in the 
asymptotic radial range to negatively curved FLRW (panel (a)) and Einstein de 
Sitter (panel (b)) backgrounds. Both models have the same central density and 
comparable density gradients in the elliptic core (the shaded area). The curves 
correspond to the initial value functions: (47r/3)pqo = 0-1 + 9-9/(1 + r^), tCq(, = 
-0.75 - 1.55/(1 + r3) (panel (a)) and (87r/3)p,o = 1-0 + 9.0/(1 + r^), Kqo = 
(1 — r-^)/{l + r^) (panel (b)). Notice that q tends as r — > oo to positive values 
in both cases (deceleration) but becomes negative (q —0.1) in the intermediate 
radial range of the model converging to negatively curved FLRW. To obtain the 
same negative value of q for the model asymptotic to Einstein de Sitter its central 
density must be about 100 times larger. 




Figure 7. Effective acceleration in the intermediate radial range: 
models converging to vacuum LTB. The figure depicts the deceleration 
parameter q (thick curve), together with fiC) and O'™), for an initial slice 
t = ti = for two models radially converging to a vacuum LTB state: the 
elliptic model asymptotic to Schwarzschild (panel (a)) and the VD hyperbolic 
model of class (A) of section 8 (panel(b)). The parameters of the models are: 
(47r/3)p,0 = 10/(1 + r3), ICqo = 1-8/(1 + r*) (panel (a)) and (47r/3)p,o = 
10/(1 + r^), K.qo = -1.25/(1 + r5/4) (panel (b)). Notice that q is negative in 
the intermediate and asymptotic radial range, and in the limit r —> oo we have: 
q — > for the elliptic model and q — > < 0, with given by ( |54[ |, for the 
hyperbolic model. 
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(i) Elliptic and hyperbolic models converging to Einstein de Sitter. These mod- 
els must be MD in the classification ( 45 ) and are characterized by the asymp- 
totic limit riq^ — 1. Conditions for a mid-range q < are far more restrictive 
for elliptic models because f^g > 1 must hold in the whole radial range, and 
thus n'^'") ^ 1 can only occur if 1 -|- (5^''^ (very large density gradient: 
scenario (38a)). The conditions are similar but less stringent for hyperbolic 
models because < fig < 1 holds in the whole radial range. 

(ii) Hyperbolic models converging to negatively curved FLRW. These must be 
G models in the classification ( 45 ) . The conditions for a mid-range effective 
acceleration are far more favorable than for models asymptotic to Einstein 



de Sitter. Since < l^o 



< 1 holds, then it is easier to obtain ilg ^ 1 in the 



intermediate radial range if we choose fig^ ^1). 
The comparison between the deceleration parameter q in hyperbolic models 
with similar density profiles that are asymptotic to these FLRW backgrounds 
is displayed in figures 6a and 6b. See sections 13 and 15 for further discussion on 
the mid-range effective acceleration in models asymptotic to FLRW. 

Models converging radially to LTB vacuum. The conditions for a mid-range 
effective acceleration are easier to meet because of the asymptotic low density: 

(i) Elliptic models. Since we have fig > 1 for all r, then OS™-'^ <^ 1 can only be 
achieved if 1-f J^''^ ^ 1 holds (scenario (38a)). This is far more likely to occur 
(from (13), (15) and ( [l7| ) in the models of case (D)(ii) in section 8, in which 
Pg — >■ sharply decays as (a = 3), so that: S^f'^ —J- —1, M — M^^,,„ and 
p <^ Pq hold (convergence to Schwarzschild |S7j)- Figure 7a shows an effective 
acceleration in the mid radial range of a complete slice of this elliptic model. 
As shown by the figure, we have q —0.1 with q becoming less negative as r 
grows so that q asymptotically (in agreement with the results of section 
9). 

(ii) Hyperbolic models. The fact that we have < given by (54) and (57 1 in 
the VD models and in some of the G models with F ^ oo clearly indicates 
that q < should hold in this range and tend asymptotically to (54) and 
(|57l 



As depicted in figure 7b, there is effective acceleration in a complete 
slice in the intermediate radial range of a VD model corresponding to the 
case (A) of section 9, with q « —0.1 and q becoming less negative as r — )■ oo, 
reaching the asymptotic value q ~ —0.03 obtained in (53) and (54 1. For 
other models the possibilities for q < depend on their being MD, CTor VD 
in the classification given by ( 45 ) . In general, since rjg ^ 1 and 1 -h 5'^P'> > 
hold asymptotically for MD models, the conditions for them are analogous 
to those of models converging to an Einstein de Sitter background. 



11. Effective acceleration near the expanding and collapsing singularities. 



We examine separately the cases where these singularities are simultaneous and non- 
simultaneous. 
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Figure 8. Expanding and collapsing singularities. The figure depicts the 
curves [thb{'i'),r\ and [tcoii{?'), r] in the {t,r) plane, which provide the coordinate 
representation of non-simultaneous expanding (big bang) and collapsing (big 
crunch) singularities for regular hyperbolic (panel(a)), elliptic (panel (b)) and 
mixed elliptic/hyperbolic (panel (c)) models. The shaded regions correspond to 
values of t allowing for complete slices '^7~[t], given hy t > tbb(O) for hyperbolic 
models and tbb(O) < t < tcoii(O) for elliptic and mixed models. The letters "C" 
and T" denote individual complete and incomplete slices, so that the latter are 
only defined for r > Ts- For the elliptic and mixed elliptic/hyperbolic models 
the locus of the bounce surface where Hq changes sign is depicted by the curve 
imax(''); whlch ixi general (see [57]) diverges as r — > oo (elliptic models) and as 
r — >■ r;, (mixed models). The comoving radius r = rj, in panel (c) marks the 
boundary of the elliptic core, so that Kq > for < r < rj, and Kq < for 
r > rf, (hyperbolic exterior). Notice that in this case, tcoii is only defined in the 
elliptic core, so that all layers collapse in this region but expand perpetually in 
the hyperbolic exterior. 



11.1. Effective deceleration near a simultaneous big-bang. 



Regular hyperbolic and elliptic models exist in which ^bb = holds with tg constant 
and tbb given by (A. 15 1. In elliptic models we can have either tl^^ = or a simultaneous 
collapse i'^ji = (but not both) [^nilMlllSISHllH]- The behavior of scalars approaching 
these singularities can be examined in the limit t ^ tg for a fixed arbitrary r. The 
following asymptotic limits hold in this limit when = (the same limits hold for a 
simultaneous collapse, see [SHUHI] for detail): 

1, S'-P^ 0, ^ 0, (70) 

while imposing t[^^ = in ( |A.15 ) and expanding (A.6), (A.8|, (A.9I (expanding phase) 
and (A. 12 1 for a w (which is equivalent to i w tg) we obtain 

^ (71) 



r(p) 



'S{t - tg 



where 2mqQ — (87r/3)/9go- Since (70) and ([tT] ) imply V. « Hq, then considering (A.l) 
and (A. 3) the integrals in the definition (ItI) of Hp are 

2m,o(l + 4''^)r'df 



■ y.R'^R'dr 
R^R'df _ 



na'^VrMf 



3 



(t-tg) 



9 



(t^tgY^ 



J" 

2m,o(l + 4^V2df 



(72 a) 
(726) 
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Hq for t : 



ts- Hence, we have Hp/Hq — > 1 as — > tg. As 
fjC^' ^ and ^ 1/2, and 



and clearly imply that Hp 

a consequence, and considering (70), we have in (311 
thus q — 1/2 > holds in this limit. Therefore, there is effective deceleration near 
a simultaneous bang singularity in hyperbolic and elliptic models. The same result 
occurs near a simultaneous crunch (t'^jj = 0) in elliptic models with this feature. 



11.2. Non-simultaneous singularities. 

As shown in figure 8, the coordinate locii of non-simultaneous curvature singularities 
are curves [ibb (bang) and [taau{r),r] (cr unch) in the {t, r) plane, with the bang 
and crunch functions, and ioon, given by (A. 15) and (A. 16). The slices ^T[t] of 



hyperbolic models marked by i < tbt,(0), and in elliptic or mixed models those marked 
by i < tbb(O) and t > tcoii(O) (shces "I" in figure 8), do not cover the radial range 
< r < Ts, where Ts is given generically by t = ibb(''s) or t = icoii('^s) for these values 
of t. The integrals in the definitions ([7| and ^ for these slices must be evaluated 
in the reduced domains 'dsi'''] = {r | r<; < r < r} C 'i?[r]. Notice that even if the 
scalars {p, %, /C} diverge as r — r^, their corresponding integrands in ([7]) and (|8]) are 



bounded (see (75) further ahead) 

Since we have < — ibb ~ ibb(''s)(^ 
with r « Ts, the limits t — >■ ttb and t — 



- rs) and t — t^^u ~ ^'oii(^s)('' ~ ^s) in domains 
fcoii are (in general) equivalent to limits r — >■ Ts . 
However, this is not the case for domains in which tl^^{rs) ~ and i'„ii(rs) w 0, such 
as very near the center [rg w 0) and in the radial asymptotic range (r^ — ^ oo) of 
LTB models asymptotically converging to FLRW or Milne (see ^7\). In these cases 
the conditions of a simultaneous singularity that were examined previously apply, and 
thus we have in them effective deceleration even close to singularities. 

Excluding domains in which tl^{rs) « and t'^^^^{rs) « hold, the following strict 
limits hold as r —> rs in both t — >■ ttb and t — >■ icon (see [58]): 



lim S'^P^ = -1, 



lim Qq — 1. 



(73) 



bmce imphes ^ and q -> n^'^'> in tu, we obtain a nonzero limit of q 

as r — >■ Ts only if the limit of the ratio 1 — Hp/T-Lq is nonzero, as both "Hp and T-Lg 
diverge positively (negatively) as t ibb {t icon)- Using the definitions ([7| and ([8| 
for A = 'H in the integration range ?9s[r] together with (A.3) we obtain the following 
rigorous general result: 

/; R^R'df /; UF-^R^Rdr 



lim 



nq{r) 



lim 



1 



lim q = 0, 



(74) 



were we applied L'Hopital's rule to each quotient above (since each one of the integrals 
vanishes separately as r — )■ r^). 



The strict limit ( 74 ) provides no information on the sign of q for domains in slices 
'^T[t] that intersect the singularities for which r > rg (marked by "I" in figures 8a-8c). 
For domains that are 'close' to the singularities in which r « we can evaluate the 



sign of q by using the analytic solutions (A.8)-(A.9) and the scaling laws (A.3), (A. 5) 



and (A. 12) to obtain the following expansions for the integrands in the definition ([7| 



This result was inferred but not proved in reference 1521 . 



-KdrsWir-r, 



(75) 
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which are bounded as r —> rs (the minus sign is consistent with the regularity 
conditions (A.13)-(A.14)). Notice that for r in the case t k, t^^^ we obtain 

the same expansions (75 1, with i'^jX^s) instead of ^(.^(^'s) and without the minus sign 
(consistent with the regularity conditions ( A.14[ ) that require t[^^^ > 0). Inserting the 
first order expansion S ~ S{rs) + S'{rs){r — r^) for S = mqor^/jF = M/J^ into the 
integrands ( 75 ) and evaluating the integrals in Q and ([s]) the ratio Hp/'Hq becomes 
at first order in r — r » 



T~La 



6J" 



(76) 



where , denotes evaluation at r 



1 + in \k.b\ with (5, 
in (31 1 at first order in 



(p) 



c(p)i 



5. Considering that « 1 + 0{L) and expanding 
together with ( 76 ) , yields a deceleration factor q 



3A/' 
2M 



rs) 



2[1-/C, 



90» 



9(1 + 4^^) 
2r 



(r - rs), 



where we used (A.l), (A.4) and (17). Since (A. 13) and (A. 14) imply 1 



(77) 

- > 0, then 

the necessary (not sufficient) condition for q < above is just the necessary condition 
( 36 1 evaluated at r = . Considering the allowed form of and the results of lemmas 
1 and 2 (see section 5 and [M]), we have the following possibilities: 

Domains in hyperbolic models. For domains r w with very large in models 
converging to a vacuum LTB state, the asymptotic convergence forms (46)-(47) 



(see table 2) transforms (77) into the asymptotic form: 



1 



Akolr^' 



3[1 - a/3] 



ir 



0</3<2, 0<a<3, 



(78) 

Considering the values of the parameters a, (3 in table 2, the following cases arise 
for domains r « in the asymptotic radial range: 

(i) There is necessarily effective acceleration in all models in which — > 00 and 
J^' — 00 hold, corresponding to < /3 < 2 (see table 13 1. In particular, q 
can be negative and large if < /3 < 1. 

(ii) There is necessarily effective deceleration in all models in which T > 1 
and J-' hold (corresponding to /3 = 2), with the exception of VD models 
with a = 3, as in this case both terms inside brackets in (78) vanish. 

It is important to remark that the results above are fully consistent with the 
results of sections 9 and 10. However, a comparison with section 9 and table 
2 reveals that an effective acceleration may exist in slices t « tbb('"s) in some 
hyperbolic models in which this acceleration in latter times was either proven 
not to exist (MD models with T — > 00) or be very constrained (G models with 
J-" — >■ 00). While the asymptotic expansion (78 1 is strictly valid only for large rs, 
this expansion indicates that it is likely that a large effective acceleration may 
exist in domains with rg in the intermediate radial range, or even in complete 
slices t > with t « tbb- As we show in the numeric example of figure 10a, 
this effective acceleration does occur (even in complete slices t ~ tbb) in a mixed 
elliptic/hyperbolic model converging to a negatively curved FLRW state. 
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Domains in elliptic models. In domains contiguous to the center in wiiicli (361 



does not hold (see lemma 2) and section 5), we have necessarily q > 0. However, 
if there is a TV (turning value) of F at some r — r^^ (see lemma 2b), then q < 
may hold in more "external" domains with {rs,r} > rtv However, since T <\ 
holds in these domains, we have the same qualitative situation as in hyperbolic 
models with T — ^ > 1. This can be appreciated by writing the asymptotic 



expression analogous to (|78|): 

13 



1 



6 



-kor^ 



2-/3 



3[1 - a/3] 



(r - rs), 



/3>2, 0<a<3. 



(79) 



> 2 when p6| holds, the first term inside the brackets is negative 



grows. Also, (A. 15) and (A. 16) imply t' 



for every 



Since we have 
but tends to zero as r, 
r, hence the curve <coii('') will show a large deviation from the lines of constant 
t (slice t = t2 in figure 9a), and thus the condition 1 + 6^^^ ~ (which implies 
^ Q-) jg fg^j. jggg Yikeiy iQ hold in domains that approach or hit the collapse 
singularity than in domains that approach or hit the big bang (see the caption of 
figures 9a and 9b). Since q > holds for all models with a < 3 as Tg becomes 
large, it is not likely that q could be negative for smaller values of r^. The 
situation may improve if a = 3 (models converging to Schwarzschild, case (D)(ii) 
of section 9), so that q < may hold in some mid-range domains (see figure 
7a). Considering the arguments expressed above, the conditions for an effective 
acceleration near the expanding and collapsing singularities are (at best) similar 
to those of hyperbolic models with T — >■ J^^ , and far more restrictive than in the 
hyperbolic cases with — >■ oo. As shown in the numeric example of figure 10b, 
there is effective deceleration in domains hitting the collapse singularity. 



12. Back reaction and effective acceleration in structure formation 
scenarios. 

Useful toy models (known as "spherical collapse models") of cosmic structure in an 
expanding background follow from a mixed LTB configuration in which a collapsing 
elliptic core region {0 < r < r^) is smoothly matched {JCq passes smoothly from 
positive to negative at r = r^) to an expanding hyperbolic exterior (r > rf,). A 
configuration of this type is illustrated by figures 8c and 9a. Back-reaction have been 
studied for the spherical collapse model by Rasiinen [171 HH] (who did not assume a 
smooth transition) and Paranjape and Singh [50] . Contrary to the claim expressed in 
[47l |48] , collapsing conditions do not favor an effective acceleration (this fact has been 
noted in the more recent references [121 [SD]). However, an effective acceleration does 
arise in the intermediate radial range of the spherical collapse model in the example 
in figure (6a). 

The results of lemmas 1 and 2 apply for spherical collapse model configurations 
(see [54]). Since < 1 must hold in the elliptic core and > 1 holds in the hyperbolic 
exterior, it is evident that the elliptic region must have a TV (turning value) of at 
some f = rtv < ''fj, so that J- must increase in the outer layers of the core (see figure 
11a). The sign of Q[r] depends on whether the domain i?[r] lies entirely within the 
core or reaches into the exterior region. In the latter case the sign of Q[r] is the same 
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"(a) 

Figure 9. Local density near the singularities of the spherical collapse 
model. Panel (b) depicts the radial profile of logjo for two incomplete slices 

marked by t = -0.056 (label "i = fi") and t = 1629.01 (label "t = t2") that hit, 
respectively, the big bang tbb{r) and the big crunch tcoii('') singularities in a mixed 
elliptic/hyperbolic model like that of figure 6a. The parameters of the model are 
the same K.qo used in figure 6a, and rriqQ = (47r/3)pgo = 0.1 + 99.9/(1 -|-r^). Panel 
(a) shows each one of the slices in the (t, r) plane. Notice that the vertical axes 
in the plot of panel (a) has very different scales for positive and negative times, 
so that tbb has a negligible time variation in comparison with tcou. While we 
have rj'*"' — > as r — ra for both slices, for the slice t = t2 that hits tcoii('') 
the function Q^™-^ diverges as the slice hits the bounce surface tmax(''), whereas 
for the other slice O'™) <gC 1 holds for all the radial range. Since 
facilitates q < 0, then it is evident that conditions for an effective acceleration are 
much more likely to be met for slices t = t\ near the big bang than t = t2 near 
the collapsing singularity (a large central density f2qo(0) = r2''"'(0) = 100 was 
selected to emphasize this feature). 



as that of domains in hyperbolic models (lemma 1), while in the former case the sign 
is the same as in domains of elliptic models with a TV of T (lemma 2(b)). Even if 
there are TV's oi H, we have Q[r] > in all domains with sufficiently large r ^ rf, 
extending in the hyperbolic region. The same remarks apply if we have a parabolic 
(instead of hyperbolic) exterior [54, since in this case T = 1 must hold for r > r^, 
hence must have a TV and grow to reach F —1 inside the elliptic core, which means 
that Q[r] > holds already for domains reaching the external layers of the core. 

While the expansions ( 78 1 and ( 79 ) arc respectively valid in the elliptic core and 
the hyperbolic exterior, some important differences emerge in comparison with the 
pure hyperbolic and elliptic models: the big bang occurs in both regions, but the 
collapse singularity must occur only in the elliptic core. Hence, t^^-nir) only reaches 
up to r = rjj and we have icon — > oo as r — >■ rf, (compare the curve t^on in figures 8c 
and 9a with the same curve in figure 8b). As a consequence, if we consider domains 
r « with rs > rtv near the collapse singularity, we have J^{rs) < J'{ri,) — 1 and 
•7^'(^s) < ^'{fb) for all Ts, and thus the term T' jJ- in (79) cannot be arbitrarily large. 
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(a) (b) 

r 




y n 

Figure 10. Effective acceleration near the singularities of the spherical 
collapse model. Panel (a) depicts the deceleration parameter q as a function of 
r for slices near the big bang for the same model of figure 9. The complete slices 
are marked by t = —0.046, —0.4, 0, 0.1 (with to = 0). These curves extend all the 
way to r = and reach negative values only in the intermediate and asymptotic 
radial range well inside the hyperbolic exterior (r;, = 1.06). The slices hitting the 
singularity are marked by t = -0.055,-0.064,-0.096,-0.137 and reach q = 
at various r = r^. Panel (b) shows the acceleration function A(r) for the slices 
t = 1437.58, 2960.17 that hit the collapse at various r = Ts (as t = t2 in figure 9a). 
Since A(r) < 0, then q > holds for these slices and domains in them effectively 
decelerate. 



As a consequence, conditions for q < in the collapsing regime are (at best) similar 
to those of hyperbolic models with J- — )• J-^ > 1, and thus are far more restrictive 
than near the big bang. 

Another nuance in domains near the collapsing singularity comes from the fact 
that every slice hitting t = t^^^ will also hit the locus of the bouncing surface T-Lq — 0: 
the curve t = t^^^{r) in figures 8c and 9a, with t^.^^{r) — oo as r — )■ r^. As shown in 
figure 9b, q diverges at the value r — r„„^ such that t = t„„x(?'max) is reached along 
any slice ^T[t] with t > t„ax(0), which includes the slices t > t^on(O). However, this 
effect is simply an artifact of the definition of q in (|31[), as there is nothing peculiar 



in the back-reaction term in the numerator of (31 ) as Q(r) -> Q{r^„^^). In principle, it 
could be possible (if r^^^ > r^v) that q — > — oo could occur as r — >■ r^^^ in these slices, 
though this possibility is quite contrived as we also have (by construction) 17^ — >■ oo as 
r — ^ r^-^^^. Since t^„a becomes very large (the curve icoii('") bends upwards and grows 
rapidly in figures 8c and 9a) in the outer parts of the elliptic core (rtv < r < rf,), then 
in the slices that hit t^^n the domains [r] rapidly reach values of r away from the 
curve L^i{r) where the quantity 1 + (5^''^ may significantly differ from zero (its strict 
limit ( 73 1), hence the term fi^™) in ( 33 ) that contains the positive contribution to q in 
(31 1 can easily become large enough to offset the negative contribution given by il'"'-' 



in (32 1. This possibility is depicted graphically in figure 9b. 



Since we can only prove analytically the value of q in domains in which r — w 0, 
we tested numerically the presence of an effective acceleration in domains in which 
r — is not small in slices near the big bang and near the collapsing singularities 
(figures 10a and 10b), for the configuration depicted by figure 9b. Bearing in mind 
that q diverges as r — )■ r^^^ in slices hitting the collapse singularity (see figure 9b), we 
computed for this case the scalar A in ( 30 ) instead of q [notice that A > is equivalent 
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to q < 0]. As shown by figure 10a, slices near the big bang and intersecting it exhibit 
large effective acceleration (up to q w —0.5) in mid range domains extending into a 
hyperbolic exterior that is asymptotically FLRW with negative curvature = 0.3). 
On the other hand, as shown in figure 10b, domains hitting the collapse singularity 
exhibit effective deceleration, just as domains in all complete slices (not displayed) 
near the collapsing singularity {t ~ <^„u(0) with t < icoii(O)). These numerical results 
are consistent with the analytic results and inferences obtained in this section, which 
clearly show that the conditions for q < are much more relaxed for t « i^b than for 
t « tcoii ■ These results are also consistent with the numerical results in [55^ , as well as 

mm- 



13. Effective acceleration in the time asymptotic range. 



Looking at the effective acceleration in the asymptotic limit t — > oo for finite comoving 
radii r is only possible for hyperbolic models, as the evolution time range of dust 
layers in elliptic models is necessarily finite because they collapse in a "big crunch" 
singularity at finite t. The existence of an effective acceleration in this "late time" 
evolution was discussed in [46] . but as we show in this section the results of this 
reference are mistaken. 

The asymptotic radial dependence of pq, ICq along the ^T[t] only depends on the 
initial conditions (the initial value functions PgOi ^gOi see [S7]), which can be chosen 
freely (as long as they comply with regularity conditions). On the other hand, the 
asymptotic time dependence of these scalars along comoving worldlines is completely 
determined by the solutions (A.7)-(A.8| of the Friedmann-like equation (A.6) (or ^ 
in the standard variables). As a consequence, the conditions for the fulfillment of 
(31 1 in the asymptotic "late time" regime are far more restrictive than those in the 



asymptotic radial regime. 

Since we are interested in the limit t 
of hyperbolic models, we can consider the functions Xi{r), yi{r) in (|A.8[) as finite 



oo for comoving observers (r finite) 



constants as we look at the late time regime from the asymptotic expansion of (A. 8 1 



on the variable a (the scale factor defined in (A.3|). Up to the leading term, (A. 8 1 



yields the following late time (t ^ tg) forms for the scale factors a and F: 



= 1 



F = 1 



\iCqor'it~to), 

3\JCqo\^/H^O^^ jt-to) 

2a 



1 



(80) 

O(toA), (81) 



The scale factors a^ and F^ approach the exact forms of a and F associated with 
the solution (A.7) for the vacuum LTB models. Therefore, for sufficiently large t 



along comoving worldlines all non-vacuum hyperbolic models should converge to these 
models. As a consequence, applying (|80|)-([8T|) to the scaling laws (flSl), ([T9|) and (|A.5 1 



yields the following asymptotic forms (r = (t — fo)/^o ^ !)• 



T-Lq 



u 



47r ^ 



2\K.qQ 



3/2 _2 



47r 
3 



2|/Cqol 

Pgo(l + 2(5, 



(p) 



0(t-^ 

2"0 ) 



)t3 



2|/C,o|(l + 4''^) 



2|/C,oP/2(l 



{K) 



(82 a) 



(826) 



(82c) 
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It is straightforward to prove that A — as r — oo. The hmits H 
as T 



oo follow from (82&) and (82 c 



and H ^ 

and these limits imply T-Lp — > (the proof is 
analogous to that of the limit as r — > oo in (40&I, replacing r by t). Since "H, "Hg, "Hp 
tend to zero, we have Q — >■ 0, and since p — > follows from ( 82a[ ), then A — > holds. 
However, it is not trivial to estimate if conditions exist (as claimed in [3^) for A to 
be positive as it tends to zero in this limit. To address this issue we need to obtain 
the time asymptotic form of Hp, which follows from integrating (over r) the time 
asymptotic forms of the integrands of ([7| obtained from (80 1, ( [8l] ) and (82c I, leading 
to: 



1 , /;[ 



^PqO( 



which, by looking 
fluctuations H — H 



1 + 25, 



(p) 





(/C) 



)-2|/C,o|(l + 4''^)]^ 



2"0 



+ 0(t3), (83) 



the time asymptotic forms (|82&|)-(|82c|), 

7 



at 

and 7i 



Hg both decay (at least) as r 



implies that the 
so that Q(r) decays 



as T~ . The sign of Q(r) in the late time regim e will depe nd o n the values of the 
coefflcients of order in the expansions (826|-(82c) and (826) (or on coefficients 
of higher order on if the coefficients of order r"'^ vanish for a special choice of 
initial conditions). Therefore, for all initial conditions (compatible with regularity) 



and irrespective of the sign of Q(r) at late times, we have from (29) and (82a): 

p(r) « 0(t-3) » Q{r) « 0(t-4), (84) 

and thus, as long as p > holds for finite times (non-vacuum models), we have 
necessarily A < (effective deceleration) as A ^ (or, equivalently, q > as q — )■ 0) in 
the late time regime, thus violating conditions ([25|, ( pO] ) and (jsTJ in this regime. We 
remark that this result is valid even if density is very low (but nonzero), because as 
small as it may be for any finite time, it will still decay slower (as t~^) than Q(r). As a 
consequence, it is evident that the results of reference |46j are mistaken. We may have 
A > as A in the late time regime only for the vacuum LTB models (see figure 
2) or for vacuum regions (such as the vacuum central region and the Schwarzschild 
exterior of section 7.2). 



14. Back— reaction, effective acceleration and spatial curvature. 

It is evident from table 1 that the conditions for Q[r] > are far less restrictive if 
JCq < (hyperbolic models) , though a positive back-reaction is also possible if /Cg > 
in domains of elliptic models in which J-' > 0, so that and ]Cq is decreasing rapidly 
along r (see section 5). While the compatibility between Q[r] > and a positive spatial 
curvature that is decreasing has been highlighted previously [48], it is important to 
discuss the following caveats. 



The "quasi-local" spatial curvature, K-q = in (13), which distinguishes 



elliptic {ICq > 0) and hyperbolic {ICq < 0) models or regions, is a different function 
from the averaged ((/C) — JCp = ^TZp/d) and local (/C — spatial curvatures, the 

latter given by Hence the signs of the three curvatures may not always coincide 
in a given domain. The averaged spatial curvature is always positive/negative in 
domains of elliptic/hyperbolic models. This follows in elliptic models from the fact 
that regularity conditions require ICq < to hold in these models 58J, so that (from 
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Models 
converging to 


limit of J- 

as r — > OO 


type 


Parameters 

a, /3 


Effective 
acceleration 


Vacuum LTB 


^ -^=0 > 1 


MD 


/3==2, 0<a<2 


No 


G 


/3 = a = 2 


No 


VD 


/3 = 2, 2<q:<3 


restricted 
(q^O) 


— > 00 


MD 


0<^<2, a</3<2 


No 


G 


< ,3 = a < 2 


restricted 
{Vlq^ < 0.065) 


VD 


0<^<2, /3<a<3 


Yes: q ^ < 


Milne 


— > 00 


VD 


,3 = 0, 0<a<3 


restricted 
(q^O) 


open FLRW 


— > 00 


G 




No 


Einstein 
de Sitter 


—5- 00 


MD 


^ = 2, a = 


No 


0<^<2, a = 


No 



Table 2. Effective acceleration in the radial asymptotic range of 
hyperbolic models. The "type" classification "MD" (matter dominated), "G" 
(generic) and "VD" (vacuum dominated) is given in equation | |45| . The 
parameters a and /3 in l |46| -| [47] l characterize the asymptotic convergence forms of 
pqO and K-qO, with the asymptotic forms of M and given by ( |48| . All hyperbolic 
models comply with Q[r] > in their radial asymptotic range (see lemma 1). 



Models 


limit of J- 


type 


Parameters 


Effective 


converging to 


as r — 


> 00 




a, 13 


acceleration 


Vacuum LTB 


^ -^00 


< 1 


MD 


/3 = 2, 0<a<2 


No, (Q[r] < 0) 








G 


(3 = a^2 


No, (Q[r] < 0) 




^-^^ 


= 1 


MD 


(3>2, (3> a 
< a < 3 


No for a < 3, 
restricted for 








G 


2 < ^ = a < 3 


a = 3 
(q^O) 


Einstein 




< 1 


MD 


/3 = 2, a = 


No, (Q[r] < 0) 


de Sitter 


^-^00 


= 1 




/3 > 2, a = 


No 



Table 3. Effective acceleration in the radial asymptotic range of elliptic 
models. The table follows the same definitions and conventions of Table 1. All 
elliptic models listed in the table in which "Q[r] < 0" is not stated comply with 
Q[r] > in their radial asymptotic range (see lemma 2). 



(lOcI and (34) applied to /C, see also figure 1) we have IC 
hyperbolic models or regions /Cp < /C, 



> /C„ > 0, whereas for 



< 



holds when JCq tends to zero (the only 
situation in which a change of sign of /Cp could occur) . As a consequence, the sign of 
Kq determines the sign of /Cp, but not (necessarily) the sign of the local curvature /C. 
This follows from the relation: 



/C = /Co 



1 + ^c^) 



/Co 



2/3 



(85) 



where we used (17) and (A.5), and F is the scale factor defined in (A.3). Considering 
the regularity conditions (A. 13 1 and (85 1, we have JCq <Q => /C < 0, so that the local 




Figure 11. Spatial curvature in elliptic models and regions. The figure 
displays the radial profile (slice t = 0) of the local spatial curvature K, = 
(thick solid curve), together with the quasi-local and average spatial curvature 
JCq, JCp and the function for a mixed elliptic/hyperbolic model like that of 
figure 6a (panel (a)) and a pure elliptic model converging to Schwarzschild like 
that of figure 7a (panel (b)). For the model of panel (a) we used the same 
{4TT/3)pqo of figure 6a with JCqo = —0.3 + 1.55/(1 -I- r®), for panel (b) we used the 
same parameters of figure 7a. Notice that we have in both cases K, < in regions 
where q < holds, even if both K,q and ICp are positive. The scalars IC, fCq, JCp 
exhibit the same qualitative behavior in all complete slices in these models. 



spatial curvature is also negative everywhere in hyperbolic models, but the converse 
of this implication is not true: /C < may hold in elliptic models or regions in which 
/Cg > and /Cp > hold everywhere. 



The implication ]Cq > 



/C > is only valid for elliptic models in which 

(/c) 



is monotonously decreasing (from (36)), and it necessarily leads to i5g > —2/3 
holding everywhere (it follows from (|85[)). Since in this case Q[r] < holds (case 
(a) of lemma 2), then positive local spatial curvature implies effective deceleration. 
However, /C, > (and thus /Cp > 0) do not (necessarily) imply /C > for elliptic 

Ttv (case (b) of lemma 2), since (36) implies 



models in which has a TV at r 



Vtv) = -2/3 and < 



2/3 must hold for r > rtv (the regularity conditions 
(A. 14) do not place lower bounds on Sq^^), hence (36) and (85) imply that JC < JCq/ 3 
< —2/3, and thus the local curvature IC becomes in these domains 



holds for 6, 



(K) 



increasingly small and may become negative as /C^ > decreases. Considering the 
asymptotic form (47) for elliptic models, (851 becomes asymptotically 

13' 



IC 



1 - 



/3>2, 



(86) 



so that /C < holds asymptotically if /? > 3, which is consistent with (5^'^-' < — 1 
holding for r > rtv, from (17), (36) and (85). In particular, it is possible to prove that 
/C < necessarily holds in regions of an elliptic core matched to a hyperbolic exterior 
at some r = rt,. Since we have ICq{ri,) — with /Cg < for r > r),, and /Cg(rb) < 0, 



then (10c 



( IO4 , ( |AT] ) and (|AJ) yield: 



IC{n) 



3 Tin) 



(87) 
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so that local spatial curvature must be already negative in the outer regions of the 
elliptic core in which /C^ > holds. The fact that spatial curvature can be negative 
in elliptic models is depicted graphically by figures 11a and lib. 

While /C > may still hold in domains of elliptic regions or models compatible 



with Q[r] > (case (b) of lemma 2) if ICq decays sufficiently slowly (as in (47) with 
2 < ^ < 3), we must necessarily have /C < for q < to occur in elliptic models. 
This follows from the parameters of the elliptic models converging to Schwarzschild 
examined in sections 9.1 and 10 (case (D)(ii)), which are the only elliptic models 
compatible with an effective acceleration (in domains in their intermediate and 
asymptotic radial range, see figure 7a and table 3). The conditions for q < in 
these models are precisely the same conditions in which < holds, namely: a decay 
of ICq with /? > 3. This is depicted graphically by comparing figures 7a and lla-llb. 
As a consequence, q < only occurs in domains (hyperbolic or elliptic) if /C < 
holds, even in elliptic models in which the averaged and quasi-local spatial curvature 
is positive. 

As mentioned in various references [5J El [SHj, the coupling between back- 
reaction and spatial curvature implies that the scaling law for (/C)[r'] — }Cp{r) in terms 



oc a 



of the average scale factor — V(p) deviates from the FLRW scaling law K. 
where a is the FLRW scale factor. It is interesting to verify the dependence of ICp 
on a-D for LTB models that comply with q < 0. Considering two classes of models 
admitting accelerating domains in the radial asymptotic range (section 9), we have: 

• Elliptic models converging to Schwarzschild (case (D)(ii)): Since — > 1, then 
V(p) cx r^, so that a-p oc r. Hence, we have ICp cx with /? > 2. 

• Hyperbolic VD models with — > oo (case (A)). Since J" (X ri-/^/^^ with < /3 < 
2, then V(p) cx r^+P/^ and av cx r(4+/9)/6. Hence, we have \ICp\ cx r-6/3/(4+/?)^ ^j^h 
< 6/3/(4 + < 2. 

Evidently, the average spatial curvature does deviate from the FLRW scaling law. In 
the elliptic case K-p decays as a power law but faster than (consistent with (/C) 
rapidly tending to zero), while in the hyperbolic case \ICp\ decays slower than a^^. 
In both cases above (hyperbolic and elliptic) we obtain the same FLRW scaling law 
only in the case (3 = 2, which corresponds to the asymptotic limit > 1 



(limit (50)) for which q > holds (cases (B) in section 9, see tables 2 and 3), though 
back-reaction is positive (at least in the radial asymptotic range) for the hyperbolic 
models (see tables 1 and 2). 

15. Final discussion. 

We have examined in a comprehensive manner the sufficient conditions for the 



existence of a positive back-reaction term (condition (29)) and an "effective" 



acceleration (conditions (30) and (31 )) in the context of Buchert's averaging formalism 
for generic open LTB models admitting one symmetry center (see [54] for a discussion 
on "closed" models with spherical topology). We have looked at these conditions by 
evaluating the involved scalars at the boundary of arbitrary simply connected spherical 
comoving domains 'd[r]. It is important to remark that we have only examined 
sufficient conditions for the existence of an effective acceleration, hence we cannot 
rule out that such an acceleration may arise in cases not considered here. 

We reviewed the following results proven rigorously in |54j (see lemmas 1 and 
2 in section 5 and table 1): back-reaction is non-negative (S[r] > 0) in the 
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Scenario 


Type of models 


Numerical 

estimate 


Reference 


Central vacuum 


Vacuum LTB region matched to 
a hyperbolic model with 
a simultaneous bang and 
converging to open FLRW 


q« -0.03 


section 7.2 
figure 2 


Schwarzschild 


Schwarzschild exterior of a dust 


not 


section 7.2 


vacuum 


ball. Hyperbolic and elliptic 
models compatible with Q > 


estimated 


figure 3 


Central void 


Under dense non-vacuum central 
region of a hyperbolic model 


q w -0.003 


section 8 
figure 4 


Asymptotic 


Hyp VD and G converging to 


q w -0.035 


section 9 


radial range 


vacuum LTB with J- oo 




figure 5 




Hyp VD converging to 


q^ 


section 9 




vacuum LTB with F > 1 








Hyp converging to Milne 


q^ 


section 9 




Ell converging to Schwarzschild 


q^ 


section 9 


Intermediate 


Hyp converging to open FLRW 


q w -0.1 


section 10 


radial range 






figure 6a 




Hyp and Ell converging to EdS 


restrictive 


section 9 
figure 6b 




Hyp converging to vacuum LTB 


q w -0.1 


section 10 




with ^ oo 




figure 7a 




Hyp converging to vacuum Milne 


not 

estimated 


section 10 




Ell converging to Schwarzschild 
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Table 4. Scenarios with effective acceleration. Further detail and 
information on the characteristics of these scenarios are found in the appropriate 
sections and figures. "Hyp" and "EU" stand for hyperbolic and elliptic, "EdS" is 
short for Einstein de Sitter, "open" FLRW refers to negatively curved FLRW, 
"VD" and "G" correspond to the classification given by \45\ and "TV" stands 
for a turning value (see section 5). The numeric estimates are the maximal values 
found in the examples. 



radial asymptotic range in all open models, except elliptic models in which T is 
monotonous, so that spatial curvature decays asymptotically as R~^. Following 
qualitative guidelines (section 6) describing specific conditions favoring the existence 
of an effective acceleration (q < or A > 0), we proved the existence of accelerating 
domains in various scenarios in sections 7-13. These scenarios are summarized in table 
4 and the details identifying the corresponding models and regions are fully stated and 
justified in these sections. 




Figure 12. Averaged vs local Hubble rate and density. Panel (a) depicts 
the ratio 'Hp/'H as a function of r along the slice t = for the hyperbolic models 
asymptotic to negatively curved FLRW of figure 6a (NegFLRW), asymptotic to 
Einstein de Sitter of figure 6b (EdS), VD asymptotic to vacuum LTB of figure 
7a (HypVD) and the elliptic model asymptotic to Schvirarzschild of figure 7b 
(EUASchw). Panel (b) depicts the ratio logjgPp/p for the same models. Notice 
that the models HypVD and EUASchw exhibit the largest values in both ratios, 
denoting the largest deviation from FLRW conditions. On the other hand, the 
least deviation occurs for the EdS model for which Hp/H ~ 1 and pp/p ~ 1 holds 
for all r, while the NegFLRW model deviates from FLRW conditions only in the 
intermediate radial range. 



It is evident, from looking at sections 7-13 and tables 1-4, that the conditions 
for an effective acceleration are far more restrictive than the conditions for a positive 
back-reaction (section 5). While condition (36) is necessary for both Q[r] > and 
q < 0, the sufficient conditions for a positive back-reaction are far more relaxed 
than the sufficient conditions for an effective acceleration. Practically all domains 
in the asymptotic radial range of hyperbolic models and of elliptic models in which 
(36 1 holds comply with Q[7'] > 0, but, as shown in the various scenarios examined 
in sections 7-12 (compare table 1 with tables 2-4), an effective acceleration exists 
only in a restricted subset of these models. In general, the conditions that (on top of 
Q[r] > 0) are most favorable for an effective acceleration are a combination of negative 
spatial curvature (hyperbolic models), regions with low density (or vacuum) and large 
gradients of density and spatial curvature, with — >■ oo being a key factor in the 
radial asymptotic range. 

As argued in section 14, even if effective acceleration may also arise in domains 
where the averaged spatial curvature is positive (the intermediate and asymptotic 
radial range of elliptic models converging to Schwarzschild, case (D)(ii) in section 9), 
the local spatial curvature in these domains is necessarily negative (see figures 7a and 
lla-b). 

While accelerating domains exist in vacuum and non-vacuum central voids 
(sections 7.2 and 8, figures 2 and 4), a density void profile (p'^ > 0, S'^''^ > 0), as used 
in most applications of LTB models to fitting observations, is not necessary at all, as 
accelerating domains also exist when density has a clump profile (p^ < 0, 5'-''-' < 0), 
provided a low density region occurs from a sufficiently fast decay of the over-density 
in the mid-range and asymptotic range, with the growth of back-reaction mostly 
driven by fluctuations (gradients) of spatial curvature (see sections 6, 7, 9 and 10). 
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In fact, the estimated numerical values of q are more negative in the intermediate 
and asymptotic range associated with these clump profiles (figures 5-7) than in the 
"inner" void profiles of central vacuum or low density regions (figures 2 and 4). 

While all models radially converging to FLRW decelerate in their radial 
asymptotic range (section 7.1), the hyperbolic models and hyperbolic exteriors with 
this asymptotic convergence are compatible with scenarios in which there is effective 
acceleration, either in a central void (section 8, figure 4) or in the intermediate 
radial range (section 10, figure 6). The conditions for this are more favorable if the 
asymptotic FLRW state has negative curvature, and this fact is consistent with the 
existence of accelerating domains in the models examined numerically by Chuang et 
al in [52] that have this radially asymptotic convergence. 

The restrictive conditions for accelerating domains in models asymptotically 
convergent to Einstein de Sitter (consistent with the results of [151111] [SU]) have been 
interpreted as a signal that back-reaction effects should be negligible when "realistic" 
models are considered [U |1H • However, the contention that only an Einstein de 
Sitter background is "realistic" has been criticized on the grounds that it imposes the 
assumptions of the standard treatment of cosmic structure of the concordance model 
(quasi-Newtonian perturbative approach and asymptotic spatial flatness [IP, "35^ ) on 
configurations in which such assumptions may not be applicable (see [51 [51 [TO, 45, 66 ). 
In particular, from an observational point of view, asymptotic spatial flatness is not 
(necessarily) required to fit supernovae and CMB constraints under a non-perturbative 
approach to inhomogeneity (see [67l [68] )■ While the relation between averaging 
inhomogeneities and the fitting of observations is not well understood and requires 
further examination (see [HI IE UHl IMllMliSn] ) j the results of a perturbative approach 
(as considered in [13 IH]) seem to be insufficient to justify a forceful conclusion that 
non-linear effects associated with back-reaction should be negligible. 

It is still an open question whether models fitting cosmic observations are also 
compatible with accelerating domains (in the context of the present work). Bolejko 
and Andersson [FT] have looked at this issue by means of specific examples: two of their 
models admit accelerating domains (models 5b and 8) but fail to satisfy age constraints 
(these two models have shell crossings and are not asymptotic to a negatively curved 
FLRW state). However, the results of these authors are based on a reduced sample of 
models and do not rule out the possibility that fitting observations can be compatible 
with models exhibiting mid-range accelerating domains and converging to a negatively 
curved FLRW state (as described in section 10). 

Some of the results proved here directly contradict previously published results. 
In particular, our proof (section 13) that all domains in the asymptotic time range 
of non-vacuum hyperbolic models necessarily decelerate proves that the condition for 
the existence of "late times" accelerating domains (found by Paranjape and Singh 
[IB]) is mistaken. We have also shown (sections 11 and 12) that the existence of 
effective acceleration is very restrictive for domains in collapsing configurations, which 
disproves an earlier claim by Rasanen in [13 US] that gravitational collapse favors this 
acceleration (this was also proven to be mistaken in [50l [49ll53] ) . 

We have obtained numerical estimates of the value of the deceleration parameter 
q defined in (jsTJ (see table 4 and figures 2, 4, 5, 6, 7 and 10). It is important to 
clarify that the main aim of this article is to study the conditions for the existence 
of accelerating domains, not to obtain "realistic" estimates for this acceleration. 
Therefore, we emphasize that the parameters used in these numerical examples 
(described in the figure captions) were selected in order to convey a rough comparative 
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inference of the order of magnitude of typical values of the effective acceleration in 
the various scenarios that were examined. Since the value of the effective acceleration 
is given by q^jf defined by (24), not by q in (31|, the obtained values of q should be 
considered as an approximated estimation of the effective acceleration. Considering 
domains in the radial asymptotic range of the VD hyperbolic models of case (A) in 
section 9, the definitions of q^jf and q in p4| ) and (31 1 yield the ratio 

2 



q A "Hp A 



q^jj. y^.'hL^ Ap 



2 



(1+0' 



where A is defined in ( [30| ) and ^ is the constant appearing in (C.9) with n = j3/2 and 
< /3 < 2, so that 7^^/W| « 1 + 0(10^2)^ However, Ap > A (considering from (IsS that 
A' < and A' — >■ must hold as r ^ oo) and Ap and A must have the same asymptotic 
limit (from (40&)). Hence, the numeric values of q^^jf and q in the asymptotic radial 
range are for these VD models of the same order of magnitude. 

Excluding domains near a non-simultaneous big bang, where values up to 
q = —0.5 were found (figure 10a), the numerical values of q in the other scenarios 
range between —0.003 and —0.1, and thus are roughly consistent with estimates in a 
perturbative context (see 5 and references quoted therein). While these values are 
small, they imply a significant difference between local and averaged scalars, which 
is an immediate consequence of inhomogeneity. As shown by the curves plotted in 
figures 12a-b, the existence of an effective acceleration is a good indicator of how 
much (and in which radial range) the averaged and local basic scalars of a given LTB 
model deviate from each other, and thus deviate from FLRW conditions. This can be 
appreciated by comparing the averaged and the local densities and Hubble rates (the 
latter approximately given by the averaged expansion scalar {%)[r] — Hpir) and its 
local counterpart H). As shown in figures 12a and 12b displaying the ratios T-Lp/T-L 
and logj^Q Ppl p as functions of r, the largest deviation of Hp and pp with respect to % 
and p occurs precisely in the scenarios favoring q < 0. In the hyperbolic VD model 
case (A) and the elliptic model asymptotic to Schwarzschild in section 9 the averaged 
Hubble rate can be up to 75% (elliptic model) and 100 % (hyperbolic model) larger 
than the local rate, while for both models the average density decays at a much slower 
rate than local density (both decay to zero). For the hyperbolic models converging to 
a negatively curved FLRW and Einstein de Sitter backgrounds the averaged Hubble 
rate can be up to 70% and 40% larger than the local rate and the average density 
can be up to 10 and 3 times larger than the local density, with these values taking 
place precisely in the intermediate radial range where q < occurs (compare the 
curves for these models in figures 12a-b with the range where q < holds in figure 
6a). Evidently, a comparison of all configurations implies that an asymptotic radial 
convergence to an Einstein de Sitter background yields a lesser deviation from FLRW 
conditions in all domains. 

The possibility that the back-reaction contribution h^^'^ in the equation for the 
deceleration parameter (31 1 could explain dark energy by "competing" with the value 
« 0.7 of the cosmological constant in the concordance model was suggested by Kolb 
et al [11]. However, this proposal has been criticized in [3 by arguing that the 
condition Jl'^'^^ « cannot be reconciled with standard structure formation scenarios 
based on plausible initial conditions for an inhomogeneous cold dark matter source 
in a FLRW background (for further discussion on this issue, see [7D]). This criticism 
is consistent with the fact that the only scenario in which we found a sufficiently 
large effective acceleration (q w —0.5) is in domains near the big bang (figure 10a), 
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but such domains would be associated with early universe conditions that cannot 
be accommodated into structure formation scenarios in such a cosmological context 
(they imply an excessively large decaying mode that cannot be reconciled with the 
accepted early universe theories and phenomenology). Other scenarios compatible 
with q < (models of section 8 and figure 4 or those of section 10 and figure 6) are 
not (at least in principle) incompatible with standard structure formation scenarios in 
a FLRW background, but the maximal values q ss —0.1 are too small to explain dark 
energy. In fact, in all models converging to a FLRW state we have q — >■ > in the 
asymptotic radial range, and thus q^ can be associated with the positive deceleration 
factor go of the FLRW dust background. 

On the other hand, the models converging to a vacuum LTB state (see sections 9 
and 10) do not converge asymptotically to any FLRW background, and thus q^^f and q 
cannot converge to the observable of some FLRW model. The same remark applies 
to the relation between the Hubble rates % and T^p and the FLRW Hubble factor, 
Hq, though averages of observational parameters should follow from an averaging 
procedure in the the light cone [IS], not along the rest frames (the slices '^T[i]). 

In fact, models converging to a vacuum LTB state for which we found q^ ~ 
—0.035 or q^ =0 negatively (cases (A) and (D) of section 9) can be understood 
in terms of the notion of "globally stationary" spacetimes intruduced by Buchert in 
[3 [221 En] J and defined by the "global stationarity conditions" given by q^^ = 
applied to domains encompassing a whole slice ^T[t] (assumed compact). Since 
q = q^ff — 0, then the models of case (D) in section 9 can be regarded as exact 
spherically symmetric (but non-compact) realizations of such class of spacetimes, with 
models of case (A) being approximate realizations (since 0.035 ^1). As shown in 
section 14, the scaling law of average spatial curvature in terms of the average scale 
factor ap significantly deviates from the FLRW scaling law (it is a power law a^" 
with < n < 2), though it does not coincide with the scahng reported in [29] (its 
equation 15), which likely occurs because the slices ^T\t] in these LTB models are not 
compact. 

Another interpretation for models asymptotic to Minkowski (either asymptotic 
to Milne or to vacuum LTB or elliptic models asymptotic to Schwarzschild) follows 
from the notion of "finite infinity" ( "fi" ) , introduced by Ellis [7l] and considered more 
recently by Wiltshire O [8], which defines a scale in which cosmological structures 
approximate bound systems that are asymptotically flat. Under this interpretation, 
the Minkowskian asymptotic vacuum field of these models can be considered as an 
approximation to a large cosmic void surrounding a spherically symmetric structure, 
and approximately describing the conditions in an intermediate scale between the 
galactic cluster scale and a cosmological FLRW far field that emerges statistically by 
considering ensembles of similar bound structures and voids [3] (see also [9l fTOl [72]). 
The existence of an effective acceleration associated with an effective deceleration 
parameter q^^f « in the intermediate "fi" scale could provide important information 
on the dynamics of larger scales encompassing ensembles of bound structures. 

While the results obtained in this article are restricted by the spherical symmetry 
of the LTB models, it is plausible that they may apply, or at least provide important 
theoretical clues, in more realistic models that depart from this symmetry. Whether 
considering possible extensions of the present article to non-spherical models (for 
example the quasi-spherical Szekeres solutions [59)). or if we remain working within 
the framework of LTB models, it is necessary to obtain numerically the precise values 
of an effective acceleration in domains of models that best fit inhomogeneities of 
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astrophysical interest, as well as exploring the connection between this acceleration 
and the fitting of observational data. Evidently, the results of the present article 
provide an effective theoretical guideline for further work along these lines. This work 
is currently under elaboration, and will be submitted in the near future. 



Appendix A. 

Appendix A.l. An initial value formulation for LTB models. 

Considering the radial coordinate gauge [J] 

Ro = R{to,r) ^ r, (A.l) 

where the subindex o denotes evaluation at an arbitrary fiducial slice t — tg. The 
conventional LTB metric ([I]) can be parametrized as the following FLRW-likc form: 
r2 



ds^ = -dt^ + + r^ (d9^ + sin^ M^^) . (A.2) 

where we have rephrased the metric functions R and R' in ([T]) as dimensionless scale 
factors 



R 

Rn 



R 
r 



R'/R 



r 1 ^) 

Rq/Ro a 

respectively mark the central and shell crossing 
singularities. The conventional free parameters M and J- in ( 13 1 can be expressed in 
terms of initial value functions as 
An 3 



so that a{t,r) — and T{t,r) 



M 



(A.4) 



while the "bang time" is given in terms of pqo, ICqo (see (A.15l). Since M — M{r) 
and J- = jF{r), then (13l, (17) and (14) lead to the scaling laws 



ni 



3a3 



a 



r ' 

3 a 



IC, 



qO, 



(A.' 



with similar scaling laws for flq, 5'--'^^ following from ( 18 ) and ( 19 ), and the local scalars 
S ^ {p, IC, H} given by = ^^(l + <5(^)). 



Appendix A.2. Analytic solutions. 

For hyperbolic (/Cgo < 0) and elliptic (/Cqo > 0) models the scale factor a can be 
obtained by the implicit solutions of the Friedman equation (A. 6) (equivalent to &)'■ 



vacuum LTB: a = 1 + \ICqQ\^/^{t - to), (A.7) 

hyperbolic (non-vacuum): j/o {t ^ ^o) = Zh{xo a) — Zh{xo), (A. 8) 

{Ze (xq a) expanding phase 

(A.9) 
2tt — Ze {xq a) collapsing phase 

fj This gauge freedom to set the radial coordinate follows from the fact that the metric ([ij is invariant 
under an arbitrary rescaling r = r{f). This specific coordinate choice is adequate for open models 
(hyperbolic and elliptic) for which Ri can be any monotonous function such that Ri{0) = 0. 
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where xq = |/Cqo|/'7igo, 2/o = |^go|^^^/"^90 with rn^o = (47r/3)pqo, and Z/i and Ze are 
u ^ Zh{u) = M^/^ (2 + u)^/^ - arccosh(l + u), (A.IO) 
u ^ Ze{u) = arccos(l - it) - u^/^ (2 - u)^/^ . (A.ll) 



Appendix A. 3. Regularity conditions in terms of initial value functions 

The condition to avoid sheU crossing singularities (which must hold for all (i, r) for 
which a > 0) is 



r = 1 + 3((5, 



(p) 



T~La 



n 



qa 



STiq {t — tQ 



> 0, 



(A.12) 



which follows from the implicit derivation of (A.8) with respect to r and using (A.9I 



and (17) to obtain T in terms of a'. The condition (A.12) leads to the Hellaby-Lake 



conditions [551 ISD ISS] that can be given analytically in the initial value formulation 
as restrictions on the initial value functions and their fluctuations [571 ISHl EI] ■ 



tL < 0, 5^'^^ + ^ > 0, S^rf> + 1 > 



(p) 



tL 



<0, 



t' 



thh — to 



>0, 



?(p) 



Hyperbolic models or regions, (A. 13) 

QD//D QD'/D "u ■ 1 > 0, Elliptic models or regions, (A.14) 

where t'^^^^ and t'^j, follow by differentiating the bang time and collapse time functions 

Zh{xo) ^ _, Ze{xo) 
— to , 

yo 

2tt — Zeixo) 



t.. 



yo 

27r 



(A.15) 



to + 



(A.16) 



yo yo 

The coordinate surface F = marks a shell crossing singularity only in non-vacuum 
models. In the Schwarzschild-Kruskal and vacuum LTB models, it is a coordinate 
singularity made by caustics of worldlines of test observers (radial geodesies in the 
Schwarzschild-Kruskal models). The surface a = Q does mark a curvature singularity 
in the Schwarzschild-Kruskal models, but not in vacuum LTB models (since p = 0). 
In terms of initial value functions, the regu larity condit ions f or the Schwarzschild- 



Kruskal and vacuum LTB models are just ( A.13) an d (A.14) with 5^''^ = — 1 and 
S'^^ = 0. As shown in the coordinate gauge ( A.l I together with absence of shell 
crossings (A.12) is sufficient to guarantee that 

i?' > and F > 0, (A. 17) 

hold everywhere, hence a radial asymptotic range is well defined, as the radial 
coordinate is qualitatively analogous to the proper radial length along radial rays. 



Appendix B. Proof of equation (28). 

We need to prove that (yV)[r] — {W{f, r)) — for 

W(f, r) = {n{f) - (H) [r]f - {H{f) - Hp{f)f , (B.l) 

where 'Hp{f) is the p-function associated with the proper volume average {'H)[r]. 
Expanding (B.l) and applying ([t]) and ^ we obtain 



V, 



^ / [2HHp-H2]-^;df, (B.2) 
A '"J Jo 
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where we used the variance relation {{^.{r) — {H) [r]Y) [r] — {V.^ 



Inserting 



= S'H = V'p/Vp and 6p = SHp = Vp/Vp in the integrand above, and bearing in mind 
that ("H) and Hp coincide at the domain boundary x = r, leads to the desired result: 



{w)[r] = -{n)M 



Vp(r) 



[VllVp]' df ^ -mM+H'pir) = Q. 



(B.3) 



An analogous result follows for the quasi-local average acting on a scalar like W with 
( )q and "Hq instead of ( ) and "Hp. 



Appendix C. Proof of the limits (50) and (51) 



We prove these limits in full generality without making assumptions on the radial 
asymptotic convergence of the scalars in (31) (other than {F^R} oo and 
{T-Lp, T-Lq] 



0). 



and — ^ as r ^ oo 



The limit (jSOj). From ( [40&| (see also [57]) we have Up 

for models radially asymptotic to vacuum LTB. The definitions of Up and "Hq in 
([7| and ([S]) can be rewritten as 



np{r) 



nq{r) 



(C.l) 



Since J- ^ J-q as r ^ oo, then for every e > there is a value f 
integration range of (C.ll so that [J^o + < J'^^ < [J'o — e 
f > y. Using (lOdi, and applying this constraint to the definitions (|C.1[), we 



y{e) in the 
^ holds for all 



obtain the inequality 
Mr) 



nq{r)~nq{y) 



< 



^0 

As r 



nq{r)^nq{y) 



V(,)(2/) 



<np{r)-np{y) 



V(p)(y) 



(C.2) 



oo we have V{p){r) oo, V(q)(r) and J-p{r) — ^ J-^ (by Lemma 4), 
therefore, if we keep y fixed and let r oo we get V(p)(y)/V(p)(r) and 
— in this hmit, so that (C.2) becomes 

^0 



^0 



-Hqir) 



3l_ 

J"n - e' 



(C.3) 



Since e can be arbitrarily small as r grows, we obtain the result ( 50 ) when J- tends 
asymptotically to a constant jFq. Notice that the converse of (50) is not true, as 
we always obtain 'Hp/'Hq -> 1 if "H — "Ho 7^ without F tending asymptotically 
to a constant. 



The limit (51). Considering (lOd), we have 

np{r) 



nqir) 



Tp{r)J{r), 



J{r) 



(C.4) 



If — > oo, then it is straightforward to prove that J-p oo, while the limit 
J^~^ — >■ imphes that for every e > there is a value x = y(e) such that J^~^ < e. 
Hence, < J(r) < e and (as long as — > oo) we have 

J(r) < 1 and lim J(r) = 0. (C.5) 
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Since Fp > 1 for all r > and it diverges as 
that Hp{r)/nq{r 



1 can only happen if J-^ ^ and J(r) have the same rate 



of convergence to zero, which is not true in general (unless J- 
Hence 



oo, then (C.5) implies 
le rate 
finite) . 



lim Fp(r)J{r) = x 7^ 1, 

r— >C30 

but if we assume expanding models then H — )• as r — 
limit, and thus we get from (34) that Hpir) > Hq{r). 
leading to: 



(C.6) 



00 implies H' < in this 
Therefore, x > 1 holds, 



hm 1 - 



n,{r) 

which is the result given in (51). 



= i-x = -e<o 



C>o, 



(C.7) 



The asymptotic form of the constant ^ in the limit (51 ) and the asymptotic value of q 
can be esti mat e d un der the assumptions on asymptotic convergence given by the power 
law forms (46l~(47). Since we assume the generic convergence law 

H ^ Hq r~", where th e va lu es o f n and "Hq depend on the asymptotic form of and 
J(^) that follow from ([46|-(47). Evaluating the asymptotic forms of the integrals in 
the numerators of ([7| and ([8|, we obtain after some tedious algebra: 



Hgir) 3(4 + /3-2n) 
which coincides with ( [5l] ) and ( |C.7[ ), since F 



< 0, 



(C.8) 



00 in ( 48 ) with fcg < requires 



< P <2. Notice that the hmit 1 
by taking P 



Hp/'Hq — in the cases T — >■ in ( 50 ) follows 



2 (see (|48|)), while the limit 1 




y.p/'Hq — of models asymptotically 



FLRW and Milne is recovered if n 
form 



The constant ^ in the limit (51) takes the 



n{2 - P) 
3(/3 - 2n + 4) ^ 




(MD models) 
(G models) 
(VD models) 



(C.9) 



where we used the asymptotic forms (46) -(47 1. 
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